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Abstract 

The  purpose  of  this  research  Is  to  study  the  hydrology  of  watersheds 
less  than  200  square  miles  throughout  the  State  of  Indiana  and  to  propose 
methods  of  estimating  the  runoff  from  these  watersheds 0  The  study  will 
include  four  parts:  a  statistical  analysis  of  existing  data,  a  geomor- 
phological  study  of  selected  watersheds ,  a  study  of  small  watershed  hydro- 
graphs  ^  and  an  analysis  of  overland  flowQ  The  present  report  summarizes 
the  work  done  from  September,  1959>  to  date  on  the  first  two  parts* 

The  statistical  analysis  has  been  made  by  means  of  extreme  values 
methodo  The  25  ysar=in3tantaneous  peak  flow  has  been  determined  for  all 
the  gaging  stations  with  a  drainage  of  less  than  200  square  miles  for 
which  the  UoS.G.So  publishes  records,,  Among  these  watersheds^  for  all 
of  which  sufficient  topographic  data  were  available,  five  geomorphological 
parameters  were  calculated:  the  watershed  area*  the  mean  relief g   the 
main  stream  slope,  the  stream  density ,  and  the  shape  factor »  A  multiple 
correlation  was  made  between  these  variables  and  the  instantaneous  25 
year  peak  flowo  A  chart  is  given  to  obtain  graphically  the  peak  flow 
knowing  the  five  geomorphological  parameters „ 


INTRODUCTION 

in  the  design  of  hydraulic  structures  or  in  the  study  of  many  hy- 
draulic projects,  a  great  and  difficult  problem  to  the  designing  engineer 
is  to  estimate  accurately  the  amount  of  runoff  that  can  be  expected.  In 
the  ease  of  highway  drainage  structures,  the  designer  is  interested,  mainly, 
in  the  maximum  runoff  resulting  from  a  storm,  The  ideal  highway  drainage 
structure  should  have  enough  capacity  to  drain  the  flood  runoff  safely 
without  causing  any  damage  to  the  structure  itself  or  to  the  surrounding 
areas. 

The  determination  of  the  required  waterway  area  of  a  bridge  or  the 
selection  of  the  slse  of  a  culvert  are  problems  which  x-equire  an  accurate 
estimate  of  the  flood  discharge  that  will  pass  through  the  structure. 
The  determination  of  this  discharge  is  more  difficult  for  small  watersheds 
because  the  majority  of  them  are  ungaged.  There  are  very  few  gaged  small 
watersheds  on  which  to  base  an  estimate,  •  Particularly  in  the  state  of 
Indiana,  there  is  very  little  information  on  watersheds  less  than  200 
square  miles.  There  are  only  twelve  watersheds  of  less  than  100  square 
miles,  and  seventeen  watersheds  with  an  area  between  100  and  200  square 
miles  in  the  state  of  Indiana  for  which  the  U,S,G.S,  is  currently  report- 
ing  flows.  The  purpose  of  this  research  is  to  obtain  reliable  methods 
for  estimating  the  runoff  for  a  safe  and  economic  design  of  highway  drain- 
age structures,  serving  watersheds  of  less  than  200  square  miles. 

The  existing  methods  are  empirical  and  fail  to  take  into  account  the 

factors  upon  which  the  runoff  depends.  Runoff  is  the  water  flowing  on 

the  ground,  the  amount  of  which  varies  greatly  with  the  storm  character- 
istics  such  aa  the  storm  duration,  intensity,  distribution!  and  watershed 

characteristics  such  as  area,  slope,  shape,  soil  physical  properties9  and 

inel  geometry. 


Unfortunately,  due  to  the  lack  of  data,  most  of  these  variables  are 
difficult  to  evaluate.  The  designers  have  therefore  used  some  empirical 
formulas  for  estimating  runoff.  The  most  famous  is  Talbot8 3  formula, 
which  estimates  the  waterway  area  from  the  watershed  area.  The  formula 
is: 

A  -  c  M^ 
V#iere:     A  »  waterway  area  in  square  feet 
M  -  area  of  watershed  in  acres 

c  «=  a  cosff icient  varying  between  1/5  and  1  depending 
on  the  slope  and  character  of  the  watershed.  The 
selection  of  the  coefficient  depends,  among  other 
things,  on  the  experience  of  the  designer. 

Due  to  the  various  factors  that  affect  the  runoff  other  than  the 
watershed  area,  the  value  of  the  coefficient  c  cannot  be  accurately  deter- 
mined to  represent  all  the  watershed  characteristics .  Talbot6 a  formula 
is  unsatisfactory  for  a  safe  design  of  a  hydraulic  structure^  Therefore, 
it  seems  desirable  to  develop  more  accurate  methods  of  computation  that 
will  include  most  of  the  factors  that  can  affect  the  runoff. 

In  this  study  the  variables  that  affect  the  runoff  are  evaluated^ 
and  improved  methods  are  proposed  for  determining  the  runoff  from  small, 
ungaged  watersheds,  within  the  state  of  Indiana  for  highway  drainage  design. 
For  the  purpose  of  this  report,  small  watersheds  are  defined  as  those 
of  less  than  200  square  miles,  that  is  those  for  which  there  is  very  little 
information  in  the  state  of  Indiana,  and  larger  than  1,000  acres  so  that 
the  land  use  is  not  a  predominant  factor. 


II.  METHODS  OF  APPROACH 

Four  methods  of  approach  will  be  considered: 
lo  Statistical  analysis 
2»  Geomorphological  analysis 
3o  Hydrograph  study 
4,  Overland  flow  hydraulics 

In  the  statistical  analysis  runoff  records  are  examined,  and  the 
theoretical  instantaneous  25-year  flow  is  calculated  by  means  of  Gumbel's 
extreme  value  theory  for  thlrty=two  watersheds  for  which  records  are 
available.  This  is  based  on  the  idea  that  the  best  ivay  of  estimating 
design  discharges  is  from  discharge  records  themselves.  This  has  the 
advantage  of  eliminating  the  need  of  relating  rainfall,  and  runoff  by  means 
of  a  number  of  variables  difficult  to  evaluate.  The  use  of  25-year  fre- 
quency is  taken  as  a  reference  frequency  which  approximates  the  expected 
life  of  highway  drainage  structures,  (1), 

A  study  is  made  to  evaluate  the  geomorphological  factors  that  can 
bs  related  to  the  stream  flow  in  order  to  estimate  the  runoff  from  ungaged 
watersheds.  If  the  watersheds  are  not  too  far  from  one  another  and  the 
storm  patterns  are  assumed  constant  over  such  areas,  there  must  be  a 
definite  relationship  between  the  geomorphological  factors  of  watersheds 
and  the  stream  flow.  Multiple  correlations  are  used  to  find  the  relations 
between  the  geomorphological  characteristics  and  the  25=year  maximum 
instantaneous  flow.  The  assumptions  made  in  the  above  approach  are  not 
absolutely  true  because  the  storm  patterns  are  not  constant  and  uniform 
over  the  gaged  and  ungaged  watersheds.  The  error  will  depend  on  the 
deviation  of  storm  patterns  among  the  several  watersheds. 


The  hydrograph  analysis  includes  the  evaluation  of  the  storm  runoff 
hydrograph  from  the  several  watersheds  and  checking  the  assumptions  made 
in  the  geomorphologisal  study.-  The  infiltration  and  detention  losses  of 
xvater  before  the  runoff  process  occurs  will  be  evaluated 0  Unit  hydro- 
graphs  and  instantaneous  hydrographs  will  be  obtained  from  which  synthe- 
tic hydrographs  for  ungaged  watersheds  will  be  established. 

The  last  approach  is  the  study  of  overland  flow  hydraulics  which 
is  based  on  the  hydrodynamics  of  sheet  flow  with  lateral  inflow.  The 
differential  equations  for  this  kind  of  flow  resulting  from  rainfall  on 
smooth  plane  surface  were  derived  from  the  general  equations  of  motion 
and  continuity.  With  the  appropriate  initial  and  boundary  conditions, 
they  can  be  solved  mathematically  as  a  propagation  problem  by  the  methods 
of  characteristics.  The  resulting  overland  flow  hydrograph  at  any  par= 
ticular  point  on  ths  plane  can  then  be  obtained.  It  is  hoped  that  with 
certain  simplifying  assumptions  that  this  kind  of  approach  can  be  used  to 
estimate  runoff  from  a  small  watershed  f@r  highway  drainage  design. 

The  present  progress  report  is  concerned  with  the  work  done  up  to 
date  which  includes  the  first  two  approaches. 


IIIo  REVIEW  OF  LITERATURE 

The  literature  review  of  the  statistical  and  geomorphological  study 
of  runoff  determination  can  be  divided  a3  three  parts: 

(1)  Statistical  theory  of  flood  frequency  analysis-~flood  forecasting 

(2)  Geomorphological  factors  which  effect  the  amount  of  runoff 

(3)  Empirical  formulas  for  estimating  runoff —relating  to  the  basin 
characteristics 

lc  Statistical  Theory  of  Flood  Frequency  Analysis 

In  1940s  Geyer  (2)  derived  a  mathematical  expression  relating  the 
flood  magnitude  to  the  exceedance  interval „  The  expression  derived  is 
as  follows:  ^ 

X-LR'* 

Vihers:    f   •  flood  magnitude  in  sec-ft0  that  has  an  excesdance  interval 

Of  (t)o 

C  ■  exceedance-interval  in  years 
L.  tRtfc-   constants  for  any  stream 
Expressing  the  above  formula  in  the  logarithm" s  form: 

/oj  Y  -  fog  L   -   i  fo$R 
/o$L-  la3Y    ■    £*  /oj  R 

/oa  f/og  L-/oa  Yj  ■   f</o$  £  +  /oj  locj  R 

It  is  of  the  form  Y  •  mx  +  b    *  which  is  the  equation  of  a  straight  line 
when/<3^(V©^£-  ~  log  fj  ±3   plotted  against  \og  £    .     If  the  appropriate 

value  for  L.    is  assigned,  then  the  plot  of  Joq(/oqL   —  /ag  /J       versus 
log  £    should  result  in  a  straight  line0  However,  the  value  of  L>     has 
to  be  assigned  by  trial  and  error  procedure,, 


In  1941,  Gumbel  (3,  h,  5,   6,  7,  3,  9,  10)  developed  the  extreme  value 
theory  „  A  special  probability  paper  was  designed  for  flood  frequency  analysis, 
According  to  Gumbel,  a  theoretical  straight  line  can  be  plotted  on  the 
special  probability  paper  with  the  observed  data  and  can  be  used  for  future 
flood  forecastingc  This  will  be  discussed  in  detail  under  the  section^ 
"Extreme  Value  Theory"  „ 

In  1950,  Benson  (11)  made  use  of  historical  data  in  flood-frequency 
analysis o  The  study  was  made  on  the  basis  of  annual  peak  discharges, and 
data  were  plotted  on  probability  paper  based  on  Guiabel8s  extreme  value 
theory o 

In  1951 »   Chow  (12)  derived  a  general  formula  for  hydrological  fre- 
quency analysis o  It  was  presentsd  to  simplify  many  existing  methods  of 
hydrological  frequency  analysis „  The  frequency  factor  /f*  depends  upon  the 
law  of  occurrence  of  a  particular  hydrological  event  «nder  consideration, 

If 

nK  *  <T  K 

Y  -  7  +nY 

Y  *  7  +<TK 

4p  -  /  +  ~^[  K 

s     /  /    CyK 

where 

^  Y  ■  deviation  from  mean 
<T  «  standard  deviation 
Y  ■  magnitude  of  hydrological  data 
Ths*  the  general  formula  is -X-  «  /¥-Cvfc    where  &  is  related  to  recurrence 
interval  7"  n 

In  1957j  Whisler  and  Smith  (13)  derived  a  method  of  estimating  the 
frequency  of  rare  floods „  The  drawbacks  of  using  the  record  of  annual 


psak  flows  which  may  result  in  errors  of  large  magnitude  in  the  determina- 
tion of  the  recurrence  interval  of  the  rarest  floods  were  pointed  outo 
This  makes  the  use  of  statistical  methods  for  such  studies  completely 
inapplicable o  Therefore,  the  technique  of  using  actual  monthly  peak  flows 
as  the  basis  for  determining  the  distribution  of  annual  peak  flows  was 
introduced o 

In  the  same  year  Rows  et'al  (14)  made  flood  frequency  analysis  by  re- 
gional synthesis  for  determining  30=~100  years  flood  for  highway  design 
in  California,,  The  procedure  recommended  was  as  follows % 
lo  Define  a  region  of  homogeneity., 
2c  From  this  region,  select  gaging  stations  likely  to  combine  with 

unbiased  analogy  0 
3„  Construct  a  frequency  graph  for  each  station, 
4°  Determine  the  five-years  flood  for  each  stationo 
5o  For  each  station,  express  all  floods  as  multiples  of  the  5=yeara 
flood o 

Qj  *  R»  Q& 

6o  Assemble  the  several  sets  of  R  in  one  synthetic  set  and  arrange  in 

descending  order  of  magnitude 0 
7«  Construct  a  frequency  graph  for  this  synthetic  record  as  though 

it  were  observed  in  a  period  of  H  year  where  N  is  the  sua  ©f  the 

periods  of  record  of  the  combined  stations „ 
80  Determine  visually  the  best-fit  line  for  the  well-defined  part  of 

the  grapho 
9o  Transfer  the  slope  of  this  line  to  each  of  the  separate  station 

graphs 0 


a 

2-.     Geomorphological  Factors 

In  1947 j»  Langbein  et  al  (15)  indicated  that  the  river  floods  were  the 
results  of  many  causes.  One  of  the  primary  objectives  of  scientific  hydro- 
logy is  the  segregation  and  evaluation  of  the  causative  factors.  The 
climatic  factor  and  the  soil=vegetation  complex  are  variables  that  exercise 
their  principal  influence  on  the  volume  of  the  runoff 0  The  topography  of 
drainage  basins  is  a  sensibly  permanent  characteristic  which  influences 
mainly  the  concentration  or  time  distribution  of  the  discharge  from  a 
drainage  basin„  The  topographic  characteristics  of  drainage  basin  cited 
were:  area  of  basin,  drainage  density  area=distance  distribution^,  length 
of  basin,  land  slope,  area-altitude  distribution,  and  area  of  water  surfaces,, 

In  1952,  Strahler  (16)  expressed  the  area-altitude  relation  by  a 
hypsometric  analysis o  The  hypsometric  curve  relates  horizontal  cross-section 
area  of  a  drainage  basin  to  relative  elevation  above  basin  mouths  By  the 
use  of  dimensionless  parameters,  the  curves  can  be  described  and  compared 
irrespective  of  true  3cale0  The  area  under  the  curve  can  be  used  to  find 
total  land  mass  and  the  mean  relief  of  watersheds. 

In  1957»  Strahler  (17)  made  a  quantitative  analysis  of  watershed  geo= 
morphology  and  showed  that  the  linear  scale  measurements  include  length  of 
stream  channels  of  given  order,  drainage  density  and  relief,,  Surface  and 
cross-sectional  areas  of  basins  are  length  products.  Dimensionless  proper- 
ties include  stream  order  numbers,  stream  length  and  Mfureation  ratio, 
maximum  valley-side  slopes,  mean  slopes  of  watershed  surfaces,  channel 
gradients,  relief  ratio  and  hypsometric  curve  properties0 

In  1959  p  Benson  (18)  cited  the  following  basin  characteristics: 
drainage  area,  channel  slope,  land  slope,  tributary  channel  slope,  watershed 
shape  factor,  mean  elevation,  percentage  of  lakes,  swamps  and  reservoirs, 
and  drainage  density*  Th3  shape  factor  of  the  watershed  was  represented  by: 


% 

//? j  equivalent  of  basin  length  divided  by  basin  width. 
*-//?  t  basin  length  divided  by  drainage  area. 
£"&£,  j  the  summation  of  snail  subdivisions  of  the  drainage  area* 

each  multiplied  by  the  distance  of  travel  to  the  gaging  point. 

3 °     I&spirical  Formulas  for  Estimating  Runoffs-Relating  to  the  Basin 
Characteristics 

In  1946,  Kinnisoa(19)  compiled  a  list  of  all  the  flood  formulas  which 

included  drainage-basin  characteristics;  there  are  the  following: 

(a)  Formula  for  extreme  or  high  flood 

Where:      Q  =  discharge  in  cfs0 
C  ■  constant 
/^  -  area  in  sq.  mio 

n  »  o.s~  0.3 

There  are  many  such  formulas  listed  in  engineering  handbooks.  However^ 
because  of  their  limited  significance  and  questionable  value*  they  are  not 
enumerated., 

(b)  Numerous  formulas  for  floods  if  an  implied  frequency  were  proposed  $, 
such  as:  the  modified  llyer  formula s   the  Kuichling,  Craig,  and 
Surge  formula,, 

(i)  Modified  Mver  formula 

Q  -  tooo  -pJ7T~ 

YJhere:      G?  =  peak  discharge  in  cfs0 
*p  ■  Myer  percentage  rating 
M  m  drainage  area  in  sq0  miles. 


10 


(ii)     Kuichling  formula 

Q  ■     ffi/lp  *   2°  frequent  flood 

Q  *       IZ7QOO     .    ?./,         - rare  flood 

Where:  <j?  =  discharge  in  sec^ft/sq.mi. 

M  ~  drainage  area  in  sq„  miles 


(iii)     Craig  fonnula 


at 


:: 


Where:      Mr  "  peak  discharge  cfs 

L.  "  length  of  drainage  basin  in  miles 
&  ■  average  width  of  basin  in  miles 
/?  *  a  coefficient  ranging  from  less  than  one  to  more 
than,  two  depending  on  the  rainfall  and  topography 

This  formula^  published  in  1868,  was  an  early  effort  on  the  part 

of  engineers  to  introduce  the  physical  characteristic  of  drainage 

basin  into  the  flood  formula. 

(iv)  Burge  formula 

s\  -ft 

Qz  /soo   -£ps 

Where:     Q  «  peak  discharge  in  cfs 
ff  »  area  in  sq<>  miles 
L.  *  length  of  basin  in  Miles 

(c)  Foraulas  for  culverts  and  storm  drains  \ 

(i)  Burhli-Ziegler  formula 


Q~  Rtfc- 


11, 

Where:  (p  °  discharge  in  cfs 

ft  »  max„   rate  of  rainfall  in  in*/i 
JS  ■  average  ground  slope  in  l«/l°00' 
C  -  a  coefficient  varying  from  0,  75 

ending  on  the  type  of  land  being 
drained 

(ii)     Talbot  formula 

Where:  n  •  the  area  of  water  way  in  sq0  feet 

C  •  &  factor  depending  on  the  topography  and  varies 
from  1/5  to  1, 
//  ■  area  in  acres 

(iii)    Farming  formula  * 

Where:  rj  «*  area  of  water  way  in  sq.  ft. 

ff  «  area  in  acres. 

(d)    Formulas  including  rainfall 
(i)     Grunaky  formula 

£  p^S  www  £OrC3/F 

%  <*  discharge  in  cms 

V   -  the  critical  time  in  minutes  during  the  continuance 
of  a  rainstorm  for  the  area  undar  consideration 
within  which  the  rain  will  produce  the  max,  rate 
of  runoff. 

K   =  Bjax0  rainfc      one  hour 

C  ■  a  coef 


:m  M  1   b- 
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This  formula  was  published  in  1922  and  is  applicable 
to  California  streams. 

(ii)  Lilie  formula 

Where:     £?«  peak  discharge 

V  ■  standard  mean  velocity  of  river  in  flood 

& mp   y   annual  rainfall 

15 

f\   ■  a  factor  based  on  the  length  of  basin  and 

is  proportional  to  the  width  of  the  basin* 

(iii)  Pettis  formula 

Q  =CpW 

Where:     Q>*  peak  flood  in  cfs 

•^2-  a  rainfall  coefficient 

W '  -  average  width  of  basin  in  miles 

0  **  a  coefficient  of  drainage  basin  characteristics 

(iv)  Besson  formula 

Qm  «  Ptti  T&mfl 

Where*       Qvn*  max,  flood 

fyti  «  max,  precipitation  during  flood  period 
c^  »  factor  representing  Influence  of  ground  condition 
/     -  factor  representing  topograph 
rt   »  drainage  area 
X  ■  an  exponent 
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(v)  Switfer  and  Miller  formula 

Where:      (p  *  24  hr0  average  flood  cfe. 
-^  »  rainfall  fao; 
|^«=  width  of  basin  in  miles 
C    and  fl    »  empirical  constant 

(e)  Formula  including  total  runoff 

Boston  Society  of  Civil  Engineers"  Formula 

Q  *  Cf&  J^ 

Where:  Q  »  peak  f3.ow  cfo. 

PjT  »  siisat  based  on  drainage  basin  character- 

R  <-  total  flood  runoff  in  inch 
/5  in  sq0 

(f)  Foraslao  with  free  impressed 
(1)                  Jller  (1914)  fev 


Q    »    ft.:  J9<«     T) 

Wh&  rage  of  the  annual  24~hrc  floods  in  cfe, 

<jj  «  raa;«,  1  day  flood  in  cfs* 
G}f?\x;  tk  diach     n  cfs* 

A^  •  drainage  area  in  sq,  miles 

"7"*  *  number  of  years  in  period  considered 

C  *  coefficient  constant 
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(ii)     R,  E<,  Horton  formula  (1914) 

f  '  *°*'-s-(-7r  ) 

Where:      4J.  .  flood  -'equalled  or  exceeded*  in  the  average 

Interval  of  T   years 
/}  »  drainage  area  in  sq„  miles 

(iii)  Eo  Wo  Lane  (1926)  formula 

Where:     .2T  «  recurrence  interval 

f(    and  &  -  constant  for  the  station  in  question 

(iv)  Wt  P„  Creager  formula 


r/y^/fcg"^/*5 


'  ■  frequency  in  years 

C  »  coefficient  depending  on  characteristics 

(g)  Flood  formula  including  basin  characteristics 

Q  *  f(S.  M.L) 
for  asxo  floods  - 

Where:     Q   «  flood  in  cfs„ 

J5"  -  median  altitude  of  the  drainage  basin  in  feet 
above  the  outlet 

rl  =  drainage  area  in  sq„  mi„ 

L.  -  average  distance  in  miles  which  water  from  run- 
off over  the  basin  .cust  travel  to  the  outlet. 


tan  txrt\  t»J  »*  itoiA*  —JJm 
jo  mM  oj 
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In  1950,  Yule  (20)  developed  a  formula  for  the  waterway  area  of  bridges 
to  be  used  in  Indiana,  which  expressed  the  waterway  area  as  a  function  of 
the  2/3  power  of  drainage  area.  It  is: 

/)  -  Co" 

Where:  /f   «  waterway  area  in  square  feet 
<2  -  drainage  area  in  square  miles 
C   -  a  general  coefficient  which  varies  with  the 
watershed  topography 
C  -  0*3  to  0„7— flat  land 
C  -  0„7  to  1.3— rolling  land 
C-  1.3  to  2.2— hilly  land 

In  1955s  Bigwood  and  Thomas  (21)  developed  a  flood-flow  formula  for 
Connecticut 0     The  method  of  analysis  was  based  on  Gumbel's  extreme  value 
theory  to  develop  a  state~wide  flood  frequency  relation j  then,  the  effects 
due  to  the  important  physical  characteristics  of  each  drainage  basin  were 
evaluated,  The  formula  has  the  following  form: 

Where:      r\   «  design  ratio  from  the  frequency  curve 

Cerr&   •  mean  annual  flood  based  on  the  watershed 
characteristics 

In  1957,  Dooge  (22)  discussed  the  rational  method  for  determining  the 
size  of  trunk  sewers,  surface  drainage  channels,  highway  culverts,  and  bridge 
waterways^.  The  essential  feature  of  the  rational  method  is  that  the  runoff 
is  calculated  for  a  storm,  the  duration  of  which  is  equal  to  the  time  of 
concentration  of  the  watershed.  The  time  of  concentration  is  that  time 
taken  for  the  first  drops  of  rain  falling  on  the  most  remote  part  of  the 


«M  le  *uq  -j.ur.oi  i  irai  *<fj 
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watershed  to  travel  to  the  outlet.  The  rational  formula  has  the  form 

VJhere:  C?  -  design  peak  runoff  rate  in  cfs. 
/?  °  area  of  watershed  in  acres 
jT=  the  rainfall  intensity  in  inches  per  hour  for 

the  design  recurrence  interval  and  for  a  duration 
equal  to  the  time  of  concentration  of  the  watershed 
C  "  runoff  factor 
This  formula  was  derived  originally  by  Thomas  Mulvany  in  1351 
for  small  watersheds  less  than  5  square  miles. 

In  1959,  Benson  (18)  found  that  the  mainstream  slope  is  the  next  in 
importance  to  drainage  area  size  which  affects  the  runoff,  The  slope  for 
that  part  of  the  mainstream  between  85  and  10  percent  of  the  total  distance 
above  the  gaging  point  provides  the  best  correlation  with  flood  magnitudes „ 
The  following  empirical  formula  was  found  for  the  New  England  region: 

■Where:  Cf  ■  peak  discharge  in  cfs. 

ft  ■  drainage  area  in  square  miles 
O  ■  85~10£  slope  of  mainstream 
&,£>.€.  *  regression  coefficient  varying  with  recurrence 
interval  of  flood 
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IVc  STATISTICAL  ANALYSIS 


Gumbel's  Extreme  Value  Statistics 

The  aim  of  a  statistical  theory  of  extreme  value  is  to  analyze 
observed  extremes  and  tc  forecast  further  extremes 0  The  flood  problem 
is  a  problem  of  extreme „  However,  no  matter  how  big  a  flood  gets,  there 
will  always  be  a  bigger  one  to  come« 

A   Extreme  value  theory 

1,  The  return  period  and  the  "expected"  extremes  value 

Let       f(x)  -  initial  probability  density  function  of  original 
data 


X 

F(x)  -  cuasulatlve  frequency  distribution  function,  that 
is,  the  probability  that  the  value  of  the  variate 
is  less  than  x 
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n  *  number  of  observations 
Un  -  a  certain  expected  largest  value 
Since  l-F(x)  is  the  probability  that  the  value  of  the  variate  is  equal 
or  larger  than  x,  then  its  reciprocal 

/ 


Too   =  /-<*?*>  4-(l) 

is  called  the  return  period.  It  is  the  number  of  observations  such  that, 

on  the  average,  there  is  one  observation  equalling  or  exceeding  x<,  From 

given  f(x)  and  hence  P(x),  then  knowing  n,  we  can  find  Un  such  that  the 

average  return  period  is  n,  Its  probability  F(U  )  is  given  by 

/  n 

2a     Exact  distribution  of  extremes 

Since  F(x)  is  the  probability  that  the  value  of  the  variate  is  less  than 
a  certain  x,  then  the  probability  that  n  independent  observations  all  fall 
short  of  x  evidently  is  F^x).  This  is  the  probability  f or  x  to  be  the 
largest  among  n  independent  observations 0  In  other  words s   the  probability 
-&n  (*n)    that  the  largest  value  fall  short  of  ^   is 
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and  the  derivative      ,-/  /„  \ 

=  nF"}'^fM  ---MO 

where j/£x„)  •  F  '  f*»)   the  initial  frequency  distribution  is  the  distribution 
of  the  largest  among  n  independent  observations. 

This  formula  constitutes  the  basis  of  the  whole  exact  theory  of 
extreme  values;  all  conclusions  are  ultimately  based  on  it„  If  the  initial 
distribution  f (x)  is  known,  then  the  distribution  of  the  extreme  values 
can  be  calculated. 

3-,  Extremal  parameters 

By  recalling  the  expected  extremes  U  defined  as  equation  4~(2) 
and  introducing^  which  is  defined  by  the  equation 

£n   -  If /fain)  4~<5) 


we  can  write 


FM    '-    i-n 


-k ■  /-  Ffu«) 


n 


G*n  - 


'-R«») 


I  -  /&„)  or 


The  reason  for  introducing  the  quantity,^   is  to  see  how  the  expected 
largest  value  Un  increases  with  n..  Differentiation  of  equation  4-(2)  with 
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respect  to  n  leads  to 

Since   c/n   -    rtf&n)  or     /?  =     —7=" 

de&«        ^  /_ 

Thus  ^   measures  the  increase  of  the  expected  largest  value  with  the 
logarithm  of  sample  size* 

4   Asymptotic  solution 

Considering  varieties  which  are  unlimited,  the  probability  function 
F(x)  converges  with  increasing  x  toward  unity   at  least  as  quickly  as 
an  exponential  type«  In  order  to  obtain  the  asymptotic  distribution  of 
the  largest  value  for  the  exponential  type,  we  expand  the  initial  probability 
function  about  the  expected  largest  value  U  by  Taylor's  expansion 


4-(7) 


•   •    o 


Since       Ft**).    *   J'lf 


F«>  *  1-4,0-%  6^f*/%,)  Mi) 
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By  using  L*Hopitalfs  rule  for  large  positive  values  of  the  variate  x, 
the  probability  density  of  f(x)  and  probabilities  1-  F(x)  are  small; 
the  derivative  f«(x)  are  also  small  and  negative c  Then,  since  f(x)  - 
P«(x) 

become  indeterminate  for  largo  values 

Then,   _j6&-  *  -  J&SL  4-(9) 

will  be  approximately  valid  for  very  large  values  of  the  variatex  since 
x  -  U„c  The  equation  4-(9)  can  be  written 


n 


4-C10) 


y?  zJ. 

from  equation  4-(5),  and  equation  4-(2) 


using  the  relation  in  equation  4-(10) 


a,  -  -  Jgf  u»i 

Multiplying  the  last  two  equations: 

/?/2i;  -  -  c2rt^  4-(i3) 

Taking  the  derivative  of  equation  4-(10)>  we  have 

tw       Cf-^fJ 

by  4-(ll)  and  4 -(12). 
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ff/?#,)   »  <3 


3 

n 


4-(H) 


The  relations  4-(13)  and  4=(14)  may  be  generalized; 
if 


n/^)  w  (-tf^*  4=d5) 


holds  for  a  certain  v,  it  also  holds  for  v  +  1 
For  the  proof,  we  start  from 


Per  the  next  derivative,  we  have 

■fa)        Q-FkM "*«,)  D-Xr.il ""  &,) 

which  proves  the  statement f  If  we  write  v  -  1  instead  of  v  *  1  and  intro= 
duce 


Into  equation  4~(8),  the  expansion  of  the  initial  probability  F(x) 


4-U7) 


*.i-*t,'t<*g&<r*0 


We  introduce  the  reduced  largest  value 


rVt> 


I- -he"  4-(i9) 


This  approximation  holds  for  any  initial  probability  function  of  the 
exponential  type  for  large  x» 
The  initial  distribution  f(x) 

ffr)     *  7f  &*€:'*  4-(20) 

equations  4°(19)  and  4~(20)  lead  immediately  to  the  asyasptotic  probability 0 


-  0-4) 

a.    n  — -  «» 


J*?x;  *  <s'  -  G*f>(~e~ )  4-(2i) 


cix 

e-e?-e-^ 

A  «-*""' 

4~(22) 


The  probability  .^fx)  and  tne  distribution  0tfO  as  function  of  yc  The 
inverse  function  Y»   -  fn(-  &ts  $0(.°))  and  the  distribution 

0^;  ^  ~£ri.)0n'JfaO     as  a  function  of  iF^x)  have  been  tabulated  by 
the  National  Bureau  of  Standards „  (9)  These  tables  are  indispensable 
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for  all  uses  of  extreme  values. 

5*  Probability  paper 

The  probability  paper  is  used  to  forecast  the  most  probable  largest 

values  to  be  obtained  within  a  certain  time  and  gives  limits  within  which 

these  values  may  be  expected  to  lis  with  a  certain  probability,  It  is 

based  on  extreme  value  theory  and  the  derived  formulas: 

A  probability  paper  designed  for  extreme  value  was  proposed  by  Power., 
The  observed  variate  x  is  traced  Gn  the  ordinate,  the  reduced  largest  value 
(y)  is  traced  on  the  abscissaP  both  in  linear  scales,  the  value  of  the 
probabilities  J§?jo  -^<?"e  )  and  the  return  period  7"  *(/-|kJ 

are  laid  off  on  auxiliary  horizontal  seales0  These,  of  course ,  are  not 
linear  scales. 

From  the  equation  relating  the  observed  variate  x  to  the  reduced 
variate  y,  /(-<<„  <f^"7?  j  it  appears  that  there  exists  a  linear  relation- 


ship between  x  and  y„  Therefore,  there re tically,  the  plot  of  x  vsc  y 
should  be  a  straight  line  on  probabilit 
of  the  line  and  d/n    is  the  intercept , 


should  be  a  straight  line  on  probability  paper,  where  ^£"JJ   is  the  slope 
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B»  Application  of  Guabel's  extreme  value  theory  (9,  10) 

Gumbel's  extreme  value  theory  shows  clearly  that  there  exists  a 
linear  relationship  between  the  reduced  largest  value  (y)  and  the  observed 
extremes  (x)  on  the  special  probability  paper.  By  extending  this  linear 
relationship,  it  is  possible  to  predict  the  future  flood  with  any  desired 
frequency,,  This  will  hold  as  long  as  the  runoff  is  kept  below  that  re- 
sulting from  the  "maximum  probable  storm"  ae  determined  from  hydrometeo- 
rological  studies, 

In  applying  Guabelf3  extreme  value  theory  to  flood  frequency  analysis, 
the  following  steps  are  followed; 

lo  Collecting  the  past  flood  data 

To  collect  the  past  flood  data  from  a  watershed  and  select  the  instan- 
taneous annual  peak  flood  as  the  extreme  value  of  each  year^  list  them  in 
increasing  order  according  to  the  magnitude  of  each  flood;  the  smallest 
flood  will  be  ranked  as  1  and  the  largest  one  will  be  ranked  as  n,  which  is 
the  total  number  of  floods r 

2<     Plotting  position 

The  plotting  position  on  the  probability  paper  is  defined  by  the  two 
coordinates:  the  discharge  x  and  the  probability  of  occurrence  <f   of  a 
raaximufflo 

Gumbel's  formula  for  plotting  position  is  given  as  follows: 

jf  =  iff/  M23) 

where     -1?  -  probability  to  be  maximum 
n    »   number  of  observations 

tYi   -  rank  of  observation  in  increasing  order  (  m  »  1, 
the  smallest  one  ) 
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This  seems  to  be  a  suitable  and  convenient  formula  for  plotting 
position,  because: 

i«  All  observations  can  be  plotted, 

ii,.  It  has  an  intuitive  meaning  and  does  not  involve  complicated 
calculations o 
iii„  It  is  analytically  simple  and  general, 
iVo  The  return  period  of  the  largest  value  is  simply  n  +  1  instead 
of  n;  the  difference  is  insignificant <> 

3o  Calculating  the  reduced  value  Y 

From  Gumbel's  extreme  value,  it  can  be  shown  that  the  relation  between 
and  Y  is 

Since  Y  is  plotted  on  a  linear  scale,  it  is  convenient  for  plotting  a 
straight  line  fitting  by  the  least  square  method „ 

A  complete  table  listing  <j>   and  its  corresponding  Y  was  given  in 
Table  2,  p,  19  of  "Probability  Tables  for  the  Analysis  of  Extreme-Value 
Data*"  (9) 

4-  Pitting  straight  line 

After  having  all  n  observations  plotted  on  the  probability  paper, 

the  next  step  is  to  fit  the  best  straight  line  through  the  plotted  points » 

The  probability  paper  is  so  constructed  that  x  and  Y  are  both  linear  scales., 

The  statistical  method  ol  least  square  may  be  applied  to  get  a  theoretically 

fitting  straight  line  with  the  least  deviations q 

Since  the  reduced  value  Y  and  the  observations  x  have  the  relation 
Y  m3(y>-u)*   which  was  shown  in  Gumbel's  extreme  value  theory,  it  can  be 

rewritten  as  X  »  U*-~  /  *  which  is  the  equation  of  straight  line,  and  the  u 
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and^  are  two  constants  which  cen  be  determined  by  the  least  square  method: 

•*  "     £yz  -fey)* 
6?  ~    x  -  &:  y 


A  more  simple  way  introduced  by  Cumbel  is  as  follows: 

x*    =  ts  ¥-  ^  Y;  i  '    t „  z,  5.  *?,  s-,    •    •    •  -  .  n 


4-(24) 


-^         z* .  cr2 


where  X  n 


X    -    CO  ■*■   ^TY 


n 


2/~\2 


x*  -a*+*v£r  +te)  (y) 


4425) 


equation  4-(2ft)  -  equation  4-(25)         X*  -  (if  *jfezr)  C  Y2  -(/)  J 


we  have         Sc~    "    .  /     p    -/y)*"  4=(26) 


and 


^  *   *    ~    <^    ^  4-(27) 


5o  Prediction  of  future  flood 

The  future  flood  with  a  certain  frequency  can  be  predicted  simply 
by  extending  the  straight  line  on  the  probability  paper.  From  the  intersex 
tion  point  of  the  extended  line  and  the  line  of  equal  recurrence  interval^ 
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the  predicted  flood  can  be  read  from  the  ordinate  of  the  probability 
paper,,  An  example  is  given  •with  complete  procedures  to  show  the  applica- 
tion of  Gumbel's  extreme  value  in  flood  frequency  analysis » 


Example: 

'Watershed  No0  S»l 

Caging  Station  -  Bean  Blossom  Creek  at  Dolan,  Indiana 


Year 

Instantaneous  Annual 

Peak  Discharge  Q 

cfs„ 

Rank 

Q 

CfSo 

£ 

y 

1946 

1830 

1 

361 

0,0665 

-  0,99716 

1947 

9420 

2 

1320 

0,133 

-  0o 70181 

1948 

2110 

3 

1390 

0o200 

-  0  47588 

1949 

9060 

4 

1740 

Oo266 

-  0o  28085 

1950 

8740 

5 

1830 

0,333 

-  0,09496 

1951 

3700 

6 

2110 

0o400 

*  0.08742 

1952 

5100 

7 

3040 

0„466 

+  0.26975 

1953 

1320 

8 

3700 

0.532 

0,46027 

1954 

361 

9 

4270 

0,600 

0,67173 

1955 

1390 

10 

5100 

0,665 

0o89657 

19& 

1740 

11 

5480 

0,733 

lol6922 

195-7 

4270 

12 

8740 

0„800 

1,49994 

1958 

3040 

13 

9060 

0.866 

1.93884 

J  vm 

5480 

14 

9420 

0.932 

2o65324 

The  discharge  Q  and  the  reduced  value  y  were  then  plotted  in  pro- 
bability paper  as  shown  in  page  (39)^  A  straight  line  was  obtained  by  the 
least  square  method <, 
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C,  Flood  frequency  analysis  for  small  watersheds  in  Indiana 
1,  The  location  of  watersheds 

Thirty~two  watersheds  distributed  throughout  the  whole  state  were 
selected  to  Btudy  for  flood  frequency  analysis.  Fig,  (1)  is  a  map  of 
Indiana  showing  the  thirty-two  watersheds,  and  Table  (1)  lists  the  names 
of  the  watersheds  and  their  assigned  number o  For  conveniences,  these 
numbers  will  be  used  in  this  report  instead  of  the  name  of  the  watersheds , 

The  purpose  of  this  study  is  to  determine  the  runoff  from  small 
watersheds „  Unfortunately,  there  are  only  three  watersheds,  the  area  of 
which  is  less  than  $0  square  miles,  12  watersheds  under  100  square  miles , 
and  29  watersheds  under  200  square  miles  for  which  flow  records  are  avail- 
able from  the  U»S,G„S.  In  order  to  have  a  good  distribution  of  the  water- 
sheds over  the  whole  state,  three  additional  larger  watersheds  were  included 
for  the  frequency  analysis,  A  bar  diagram  was  plotted  in  Fig„  (2)  to  show 
the  time  period  of  records  for  each  watershed,, 

2->  The  instantaneous  peak  flood 

Instantaneous  annual  peak  flow  was  used  for  the  extreme  value  analysis 
instead  of  the  maximum  mean  daily  flow.  The  reason  for  this  is  that  the 
hydrographs  for  small  watersheds  have  much  sharper  peaks  and  the  instantaneous 
peak  is  the  most  critical  discharge  for  hydraulic  structures,,  Data  for 
this  instantaneous  annual  peak  were  obtained  from  the  U»SoGoSo  office  in 
Indianapolis,  Indiana „  The  data  are  analyzed  on  a  water  year  basis-=«-that 
is^  from  September  to  August, 
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Figure  1 
Location  of  Watersheds  for  Flood  Records  Studied 


LOCATION      OF     WATER- 
SHEDS   FOR    FLOOD 
RECORDS    STUDIED 

Figure  I 
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Table  I  -  List  of  Watersheds,  Their  Area  and  Assignee;  Number 


Watershed  ««*„«.*.«,» 

No..                    Gaging  Station  .  Watershed 

<■"— — — —=——___ -r^n,  , .--  Area  (gq  mi) 

S=l     Bean  Blossom  Creek  at  Dolan,  Ind,  X00 

S-2     Clifty  Creek  at  Hartsville,  Ind,  88c8 

S-3     North  Fork  Vernon  Fork  near  Butlervilla,  Ind,  87,3 

S-4     Hart  Ditch  at  Munster,  Ind„  69 o  2 

S=5     Salt  Creek  near  McCool,  Ind„  73,7 

Little  Calument  River  at  Porter,  Inde  62„9 

S-7     Cedar  Creek  at  Auburn,  Ind0  93^0 

S=a     West  Creek  near  Schneiderp  Ind„  46^3 

S-9     Iroquois  River  at  Rosebud,  Indc  30o3 

S-10    Bice  Ditch  near  South  Marion,  Ind„  22.6 

S=ll    Big  Slough  Creek  near  College ville,  Ind„                84„I 

S-12    Carpenter  Creek  at  Egypt 9   Indo  48,1 

1  Tippecanoe  River  at  Oswegos  Ind,  115 

2  Missiselnewa  River  near  Bidgevill«,  Indc               130 

3  Wildcat  Creek  at  Greentown;,  tod,  162 
Cicero  Creek  near  Arcadia,  Ind„  131 

5  Fall  Creek  near  Fortville,  Ind.  172 

6  Eagle  Creek  at  Indianapolis,  tod0  I79 

7  Young  Creek  near  Edinburg,  Ind0  109 

8  Blue  River  at  Carthage,  Ind»  I87 

9  •  Sand  Creek  near  Brewers ville,  tod„  I56 

10  North  Fork  Salt  Creek  near  Belmont,  tod*  120 

11  Patoka  River  at  Jasper,  Indo  257 

12  Busseron  Creek  near  Carlisle,  Indo  228 

(continued) 
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Table  I  -  List  of  Watersheds,  Their  Area  and  Assigned  Kunber  (continued) 


Watershed  Watershed 

No...  Gaging  Station  Area  (sq  ml) 


13  East  Fork  White  Water  River  at  Richmond,  Ind,  123 

14  Silver  Creek  near  Sellersburg,  Ind„  188 

15  Big  Indiana  Creek  near  Corydon,  Ind0  129 

16  Kankakee  River  near  North  Liberty,  Indo  152 

17  Singleton  Ditch  at  Schneider p  Ind0  122 

18  Deep  River  at  Lake  George  Outlet  at  Hobart,  Ind„  125 

19  Pigeon  Greek  at  Hogback  Lake  Outlet  near  Angola^,  Indo  102 

20  Laughery  Creek  near  Farmers  Retreat,  Indo  248 
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Flgur®  2 
Period  of  Record  of  Sistantanoous  Armual  Peaks  at  Gaging  Stations 


Fig.  2      PERIOD    OF  RECORD    OF    INSTANTANEOUS 
ANNUAL    PEAKS    AT    GAGING    STATIONS. 

Watershed 
NO. 

PERIOD    OF  RECORD 

♦ 

n 

1900 
1910 

1920 

1930 

1940 

1945 

1950 

1955 
I960 

S-  1 

14 

S-2 

12 

S-3 

18 

S-4 

17 

S-5 

15 

S-6 

15 

S-7 

17 

S-8 

9 

S  -9 

II 

S-IO 

II 

s-n 

10 

S-12 

10 

1 

— 

II 

2 

13 

3 

16 

4 

5 

5 

18 

6 

, 

21 

7 

17 

8 

— 

10 

9 

12 

10 

1 

14 

II 

1 

1 

12 

12 

16 

13 

, 

10 

14 

5 

15 

17 

16 

9 

17 

II 

18 

13 

19 

14 

20 

19 

*    n=    Length    of   Record    in    Years 
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3'  Examine  the  data 

In  examining  the  data,  the  nan~made  changes ,  such  as  new  reservoirs „ 
soil  and  vatervation  works,  and  changes  of  gaging  site  will  be  carefully 
studied  because  they  can  change  the  natural  conditions  of  stream  flow  and 
affect  the  peak  discharge 0  There  are  four  kinds  of  man-made  changes  to 
be  studied  for  those  thirty-two  watersheds.  They  are  discussed  as  follows: 

i„  Changing  of  gaging  site  (23) 

There  are  four  watershed  gaging  sites  vihich  were  changed „ 


Watershed 

Sit®  Change 

Watershed 
Area 

Change 
of  Area 

Hobart 

400*  upstream 

125 

0 

Brewersville 

1^7  miles  upstream 

153 

-3 

Angola 

1,5  miles  downstream 

105 

♦3 

Greentown 

230  miles  downstream 

» 

+10 

1 

Table  II 
(The  watershed  area  changes  due  to  the  change  of  gaging  site.) 


The  above  table  shows  the  four  watersheds  whose  gaging  sites 
were  changed,  causing  the  changes  of  watershed  area.  Since  the 
gaging  station  did  not  move  far  away  up  or  downstream  from  the 
old  gaging  site  and  the  change  of  area  is  small,  its  effect  on 
the  peak  flow  can  be  ignored, 
li,  Reservoir 

A  reservoir  was  built  upstream  of  each  of  the  gaging  stations 
of  watersheds  (S-l)  and  (11).  The  U.S.G.S.  reported  that  their 
effect  on  the  peak  flow  was  insignificant.  (23) 
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iiic  Over  bank  storage 

The  backwater  caused  by  overbank  storage  may  influence  the  storage- 
discharge  relation  of  the  watershed  (l6)0  Since  the  storage-discharge 
relation  was  corrected  for  this  effect,  the  overbank  flow  will  not 
affect  the  flow  record  at  alio 
iVo  Dredging 

Some  dredging  has  taken  place  in  channels  of  watersheds  (17) and  (4), 
This  mainly  affects  the  stage-discharge  relation;  however,  there 
seem  to  be  no  serious  influence  on  the  peak  discharge  measuremento 

From  a  close  study  of  these  four  nan-made  changes ,  it  can  be  con- 
eluded  that  their  influence  on  peak  flood  is  insignificant  and  can  be 
ignored  o 

4   The  historical  flood  (11,  24) 

In  plotting  the  observed  data  on  the  probability  paper,,  a  problem 
arises  when  using  the  historical  data.  The  historical  flood  is  important 
t©  the  flood  frequency  analysis,,  It  gives  mora  information  and  makes  the 
frequency  analysis  more  reliable,  especially  in  cases  where  short  time 
records  are  available.  It  is  questionable  whether  or  not  the  probability 
curves  derived  from  short  periods  of  records,  say  less  than  25  years,  can 
be  extrapolated  to  cover  the  historical  datac 

The  plotting  position  of  historical  floods  may  be  computed  for  the 
following  two  cases: 

I,  Where  there  is  one  historical  flood ,  higher  than  any  for  the  period 
of  re cord „  In  this  case  the  recurrence  interval  of  the  historical 
flood  should  be  recorded  as  equal  to  one  plus  the  period  for  which 
it  is  the  greatest o  An  example  of  this  is  the  1913  historical 
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flood  in  Potaka  River  at  Jasper,  Indiana.  The  flood  was  16,000 
cfso  The  period  for  which  it  is  the  largest  is  47  years. 
According  to  the  formula 


3&    -  nT7 


47 


The  recurrence  interval  is 


T 


— ^-~.  ^  -5*3  Y&ARS 


ii,  The  second  case  is  where  there  is  one  historical  flood  known  to 
be  the  highest  until  a  greater  flood  occurs  during  the  period  of 
record o  An  example  of  this  is  the  1913  historical  flood  in 
Eagle  Creek  at  Indianapolis,  Indiana „  It  is  19,000  cfsc  And 
the  1957  flood  is  26,600  cfs  which  is  higher  than  the  historical 
one0  The  plotting  position  of  these  cases  are  shown  in  the 
following  Table  (3): 


Floods 

n      JF     T 

28800  cfs.  max*  (1913—1959) 
19000  cfs,  2nd  flood  (1913—1959) 

47 

47 

47/48 

46/48 

48 

24 
24 

Table  III  <=  Calculation  of  Historical  Floods 


5-  Straight  line  fitting  for  future  prediction 

According  to  Gumbel's  theory,  the  observed  esrfcremes  should  plot 
along  a  straight  line  on  the  probability  paper.  The  statistical  method 
of  least  square  is  used  to  obtain  best  fit  line  with  the  least  deviation,, 
An  important  question  arises  when  the  record  of  flood  data  is  short  and 
needs  to  be  extrapolated  to  longer  periods,.  The  historical  flood  data 
can  help  to  give  a  better  frequency-flood  relation.  But,  there  will  still 
be  a  question  as  to  the  reliability  and  recurrence  interval  of  the  histori- 
cal flood  data  itself, 
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Based  on  the  above  discussions,  it  is  questionable  whether  the  least 
square  method  for  straight line  fitting  can  be  applied,  since  the  lower 
points  in  the  probability  paper  play  a  more  important  role  than  the  upper 
larger  flow  points  in  the  true  flood-frequency  relation,,  Therefore, 
those  lower  points  should  be  weighed  more  in  the  process  of  line  fitting » 
The  line  fitting  in  this  report  is  fitted  in  two  ways: 

io  For  those  which  the  plotted  points  have  already  formed  a  straight 
line  pattern,  the  least  square  method  is  used  in  order  to  get  a 
best  fitting., 
iio  For  those  which  some  of  upper  points  are  too  far  away  the  normal 
straight  line  pattern,  then,  the  line  is  fitted  simply  by  inspec- 
tion in  order  to  give  the  lower  points  more  weight  in  the  fitted 
line. 

Results 

Following  are  the  results  from  the  frequency  analysis  for  each  se- 
lected watershed.  The  instantaneous  annual  peak  discharge  against  the  re= 
duced  value  (Y)  were  plotted  in  probability  paper,  and  a  straight  line  was 
fitted  by  the  least  square  method.   By  extending  this  line,  the  future 
floods  of  different  frequencies  were  predicted  and  shown  in  Table  IV., 
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(S~l)  Bean  Blossom  Creek  at  Dolan,  Ind,    (23) 


Location,  —Lat  39014,,30,,,  long  86°29'57M,  in  SW$  sec  2,  To  9  N„,  R.  1  ¥,,  on 

downstream  side  of  right  pier  of  highway  bridge  at  Dolan,  17c 5  miles  upstream 
from  mouth o 

Drainage  area, ==100  sq  mi, 

Gage,==Nonrecording  gage  Apr,  3,  1946,  to  Sept,  27,  1951j  recording  gage  thereafter. 
Datum  of  gage  is  576,41  ft  above  mean  sea  level,  unadjusted,, 

Stage-discharge  relations-Defined  by  current=meter  measurements.  Discharge  adjusted 
for  rate  of  change  of  stage  above  5  ft.  Only  annual  maximums  adjusted  prior 
to  installation  of  recording  gage. 

Flood  stage,  —15  ft. 

Remarks,— Flow  regulated  since  April  1953  by  Bloomington  Keservoir  (capacity, 

4,640,000,000  gallons)  7j  miles  upstream;  peak  discharges  probably  not  materially 
affected , 


Peak  Stages  and  Instantaneous  Annual  Pe3k  Discharge 


1 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

I^ay  16,  1946 

13,0 

1*830 

1953 

Mar-.  4,  1953 

11,07 

1,320 

June  2,  1947 

17,8 

9,420 

1954 

Kay  2,  1954 

5  c  45 

361 

Mar, 27,  1948 

13-5 

2,110 

1955 

Apr.13,  1955 

11,63 

1,390 

Jan,  5,  1949 

17,9 

99060 

1956 

Kay  28,  1956 

12,93 

1,740 

Jan.  4,  1950 

17.75 

8,740 

1957 

Kay  22,  1957 

15o78 

4,270 

Jan, 21,  1951 

15,50 

3,700 

1958 

June  14,  1958 

3,040 

May  24,  1952 

16.12 

5,100 

1959 

Jan,  21,  1959 

5,480 

39. 
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(S-2)  Clifty  Creek  at  Hartsville,  Ind, 


Location.— Lat  39°l6825",  long  85°42«10",  in  NW£  sec.  36,  T,  10  M.,  R.  7  E,,  at 

downstream  side  of  left  abutment  of  highway  bridge,  a  quarter  of  a  mile  north 
of  Hartsville,  and  5  miles  upstream  from  Duck  Creek0 

Drainage  area<>~~83.8  sq  mi0 

Gage0~=Nonrecording  gage  Feb,  12,  1948,  to  Sept.  23,  1952;  recording  gape  thereafter  = 
Datum  of  gage  is  677*34  ft  above  mean  sea  level,  datum  of  1929<, 

Stage-discharge  relation. ^-Defined  by  current-meter  measurements  below  6,000  cfs0 

Historical  data  „ —Flood  of  1913  on  Clifty  Creek  reached  a  stage  of  about  3  ft  higher 
than  the  KcKinley  (1897)  flood  according  to  a  report  in  the  Evening  Republican 
of  Columbus,  Ind0  dated  Mar.  25,  1913=  (The  preceding  statement  was  apparently 
for  the  Petersville  area,  about  6  miles  doxmstream  from  Hartsville) „ 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


iter 
sar 

Date 

Gage 
Height 

Discharge 
cfs 

VJater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

?13 

Mar, 

25,  1913 

25,1 

1954 

May  27,  1954 

4.17 

635 

m 

Mar; 

27,  1948 

8,48 

3,710 

1955 

July  8,  1955 

6o24 

1,760 

m 

Jan, 

5,  1949 

13  o4 

8,100 

1956 

June  22,  1956 

11,10 

5,890 

?50 

Jan., 

4,  1950 

11,8 

6,520 

1957 

July  4,  1957 

9,28 

4S270 

m 

Nov0 

20,  1950 

8,9 

3,910 

1958 

Kay    6,  1958 

2,700 

?52 

Jan, 

26p  1952 

11,5 

6,250 

1959 

Jan. 21,  1959 

11,300 

>53 

Mar-, 

4,  1953 

5.57 

1,370 
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(S-3)  North  Fork  of  Vernon  Fork  near  Butlerville,  Ind, 


Location— Lat  39°02«55",  long  85°32«40»,  in  SEfc  sec,  17,  To  7  No,  R,  9  E,,  on 
left  bank^  0,3  niile  downstream  from  Muscatatuck  State  School  dam,  l£  miles 
downstream  from  Brush  Creek,  and  2  miles  northwest  of  Butlerville, 

Drainage  area ,==37 =3  sq  mi, 

Gage,~=Nonrecording  gage  Feb*  16,  1942,  to  Aug,  18,  1942;  recording  gape  thereafter. 
Datum  of  gage  is  669 ° 40  ft  above  mean  sea  level,  datum  of  1929, 

Stage -discharge  relation, --Defined  by  current-meter  measurements „ 

Flood  stage, -==-11  fto 


Peak  Stages  and  Instaneous  Annual  Peak  Discharge 


ar 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

42 

Apr^  9,  1942 

8,94 

2,560 

1951 

Nov, 

20, 

1950 

15,98 

8,030 

« 

Mar,l6,  1943 

17  79 

9,910 

1952 

Jan, 

26, 

1952 

13,18 

5,300 

1*4 

Apr  11,  1944 

12,63 

4,780 

1953 

Mar, 

4, 

1953 

10,34 

3,260 

k5 

Mar,  6,  1945 

18,72 

10,900 

1954 

Jan, 

1, 

1954 

5=58 

840 

K6 

Feb, 13,  1946 

15,95 

8,030 

1955 

Feb, 

27, 

1955 

12o05 

4,300 

k7 

June  2,  1947 

14,30 

6,330 

1956 

May 

28, 

1956 

16,23 

8,330 

IS 

Mar, 27,  1948 

15,12 

7,130 

1957 

May 

22, 

1957 

17.04 

9,080 

k9 

Jan, 24,  1949 

18,73 

10,900 

1958 

July 

22, 

1958 

7,730 

50 

Jan  4,  1950 

17  90 

10,000 

1959 

Jan, 

21, 

1959 

26,200 

U2> « 
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(S=4)  Hart  ditch  at  Minister,  Ind« 


Location.-— Lat  41°33>40",  long  87°28»  50%  in  N  l/2  see,  20,  T.  36  N,,  R-  9  tf.,  on 
left  bank  at  city  limits  of  Munster,  a  quarter  of  a  mile  downstream  from  U  3, 
Highway  41,  and  0,4  mile  upstream  from  mouth., 


Drainage  are  a  ^^=69,2   sq  mi= 


Gage o— Recording =  Datum  of  gage  is  591 „ 21  ft  above  mean  aea  level,  datum  of  1929 , 

Stage  discharge  relatione —Defined  by  current-meter  measurements „  Dredging  operations 
assumed  tc  have  occurred  between  April  1944  and  April  1945 »   and  subsequent  filling 
have  affected  high~water  rating,  Backwater  from  Little  Calumet  River  and  possi- 
bly from  cverbank  return  affects  stage  at  gage  at  times  during  periods  of  ex- 
tremely high  floWc 

Flood  Stage o --7  ft. 

Remarks, —Hart  ditch  is  tributary  to  Little  Calumet  River,  At  this  point  low  flow 

of  Little  Calumet  River  runs  we3t  into  Calumet  Sag  Channel  or  into  Lake  Michigan 
through  Grand  Calumet  River;  floodflow  at  times  runs  east  into  channel  storage 
or  through  Burns  ditch  to  Lake  Michigan, 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


'ater 
ear 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

943 

Mar,  16 j 

1943 

6,95 

2,280 

1952 

June  14, 

1952 

4,39 

1,190 

944 

Mar,  15, 

1944 

7.23 

2,420 

1953 

Mar,   15, 

1953 

3  84 

960 

945 

May      8, 

1945 

3  73 

x,270 

1954 

Mar0   25, 

1954 

4:25 

1,110 

946 

Jan-.     6, 

1946 

2,88 

780 

1955 

Oct,  11, 

1954 

7o83 

2,600 

947 

Apr,     6, 

1947 

6,17 

2,490 

1956 

May     11, 

1956 

5,27 

1,550 

948 

May    11, 

1948 

5,60 

1,950 

1957 

July  14, 

1957 

7,60 

2,060 

949 

Febo  13, 

1949 

3,00 

850 

1958 

June  10, 

1958 

960 

950 

Dec,   22, 

1949 

4o83 

1,570 

1959 

Apr.  28, 

1959 

2,670 

951 

May  11 , 

1951 

5  01 

1,430 

U5. 
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(S-5)  Salt  Creek  near  KcCool,  Ind, 


Location.,— Lat  aa35"48n,  long  87°08»40»,  in  SE*  sec,  6,  T  36  Mo,  R-  6  W  ;  on  left 

bank  on  downstream  side  of  highway  bridge ,  50  ft  downstream  from  New  York  Central 
Railroad  bridge,  1^  miles  north  of  McCool,  and  1,5  miles  upstream  from  Little 
Calumet  River, 

Drainage  area  —78,7  sq  mi, 

— IJonrecording  gage  May  5,  1945,  to  July  24,  1955;  recording  gage  thereafter 
Datum  of  gage  is  594,10  ft  above  mean  sea  level,  datum  of  1929  (levels  by 
Indiana  Flood  Control  and  Vfater  Resources  Commission)  o 

Stage -discharge  relation. —Defined  by  current-meter  measurements  below  2,300  cfs 
' Flood  stage, =-10  ft. 


Feak  Stages  and  Instantaneous  Annual  Peak  Discharge 


ter 
ir 

Date 

Gage 
Height 

Discharge 
cfs 

Vater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

45 

June 

29, 

1945 

10.88 

990 

1953 

Mar.  16, 

1953 

816 

454 

he 

June 

13, 

1946 

11-27 

1,280 

1954 

Mar,  26, 

1954 

10.48 

910 

47 

Apr, 

5r 

1947 

11,83 

1,580 

1955 

Oct,  11, 

1954 

14-12 

3,180 

48 

May 

11, 

1948 

12,3 

1,910 

1956 

Apr,  29, 

1956 

11,26 

1,280 

49 

Feb, 

14, 

1949 

9-28 

525 

1957 

Apr,  27, 

1957 

9  81 

725 

50 

22, 

1949 

12o02 

1,700 

1958 

Nov,  15, 

1957 

456 

51 

11, 

1951 

10,78 

970 

1959 

Apr,  28, 

1959 

1,200 

52 

Nov . 

14,- 

1951 

10.63 

912 
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(S=6)  Little  Calumet  River  at  Porter,  Ind, 


Location,— Lat  41°37*18",  long  87°05'13",  in  HE  1/4  sec.  34,  T.  37  K.a  H,  6  W,,  near 
center  of  span  of  downstream  side  of  highway  bridge,  three-quarters  of  a  mile 
northwest  of  Porter,  and  4*5  miles  upstream  from  Salt  Creeko 

Drainage  area,— 62»9  sq  mi0 

Gage  --Nonrecording  gage  May  5,  1945,  to  June  25,  1952;  recording  gage  thereafter 
Datum  of  gage  is  603 =48  ft  above  mean  sea  level,  datum  of  1929 . 

Stageniischarge  relation „— Defined  by  current-meter  measurements  below  2,500  cfs 
Rating  subject  to  changes  throughout  range  of  stage-, 

Flood  stage „ ==7  ft„ 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


er 

r 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

5 

June  28^ 

1945 

9,88 

2,440 

1953 

May 

23,  1953 

6,64 

521 

6 

June  13 i 

1946 

6-99 

715 

1954 

Apr- 

26,  1954 

8-32 

1,170 

n 
1 

Apr   5, 

1947 

9-42 

2,140 

1955 

Oct. 

10,  1954 

11  66 

3,no 

8 

May  11, 

1948 

9,10 

1,960 

1956 

Apr0 

29,  1956 

8.67 

1,370 

9 

Kay  20, 

1949 

6,88 

690 

1957 

AprP 

27,  1957 

7>65 

848 

0 

Dec,  22, 

1949 

8„72 

1,720 

1958 

Feb, 

28,  1958 

490 

1 

t'z-7     11, 

1951 

8„11 

1,360 

1959 

Apr. 

28,  1959 

1,420 

2 

Nov.  14,o 

1951 

7.92 

1,060 
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(S-7)  Cedar  Creek  at  Auburn,  IncL 


Location,  —Lat  41°21«,  long  85°03S  in  SW  lA  sec  29,  T=  34  N„,  R,  13  E,  „  near  center 
of  span  on  upstream  side  of  Ninth  Street  Bridge  in  Auburn  and  2  miles  upstream 
from  Peckhart  ditch , 

Drainage  area.,— =93  3q  mi. ,  appro:dmatelyo 

Gage , ~ Nonrecording  gage  July  30,  1943,  to  Sept,  30,  1953  J  recording  gage  there- 
after Datum  of  gage  is  847ol4  ft  above  mean  sea  level  (city  of  Auburn  bench 
mark), 

Stage=discharge  relation,  ™Def ined  by  current-meter  measurements. 

Flood  stage „~8  ft. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


ter 
ar 

Date 

Gage 

! -eight 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

W 

May 

1943 

9  8 

1,470 

1952 

Mar, 

11, 

1952 

8,45 

900 

Uh 

Apr,  12, 

1944 

9-3 

1,230 

1953 

Mar, 

U> 

1953 

5,80 

471 

k5 

May  18, 

1945 

913 

1,150 

1954 

Mar, 

25, 

1954 

7o57 

707 

* 

June  13, 

1946 

8.58 

916 

1955 

Jan, 

6, 

1955 

7.61 

707 

h7 

Apr*  21  j, 

1947 

9. ,02 

983 

1956 

Apr, 

30, 

1956 

8o85 

1,050 

+8 

Feb,  28, 

1948 

8,58 

905 

1957 

Apr, 

6, 

1957 

6,89 

651 

W 

Feb,.  16, 

1949 

9,21 

995 

1958 

Dec. 

20, 

1957 

540 

50 

Apr.  5, 

1950 

9,90 

1,520 

1959 

Feb, 

14, 

1959 

890 

51 

Feb.  22, 

1951 

8,65 

940 

51. 
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(S-8)  West  Creek  near  Schneider,  Ind, 


Location.  — Lat  41°12»52",  long  87029'36,(s  in  NW  1/4  NE  1/4  sec.  19,  ?.  32  N.,  R.  9  W,, 
on  left  bank  at  downstream  side  of  county  highway  bridge,  1.2  miles  upstream 
from  Singleton  ditch  and  2  3/4  railes  northwest  of  Schnsider, 

Drainage  area,—  54 , 5  sq  mi. 

Gage, —Monrecord inn  gage  July  29,  1948,  to  Dec,  31,  1951,  and  Jan,,  1,  1954,  to  June  10^ 
1956;  recording  gage  since  June  11,  1956.  Datura  of  gage  is  627.86  ft  above 
nean  sea  level,  datum  of  1929  (levels  by  Soil  Conservation  Service). 

Stage-discharge  relation.— Defined  by  current-meter  measurements , 

Flood  stage,  .--=7  ft. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


r 

Date 

Gage 

Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

Feb, 

13,  1949 

4^58 

504 

1956 

Feb  25,  1956 

5*42 

710 

Dec 

22,  1949 

6.  56 

1,050 

1957 

July  13,  1957 

7  02 

1*250 

Feb,. 

19,  1951 

5,52 

738 

1958 

June  9,  1958 

794 

Mar- 

25,  1954 

6.10 

1,000 

1959 

Apr.  28,  1959 

1,200 

Oct. 

10,  1954 

8.06 

1,840 

53, 
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(S-9)  Iroquois  River  at  Rosebud,  Ind, 


Location. —Lat  41°02%  long  87°11' ,  in  SV.'  1/4  3ec,  24,  T«  30  N.,  R.  7  W.,  100  ft 
downstream  from  bridge  on  county  road,  half  a  mile  north  of  Rosebud ,  half  a 
mile  downstream  from  confluence  of  Swain  and  Dexter  ditches^  1.5  miles  upstream 
from  Davidson  ditch,  and  2  miles  east  of  Farr, 

Drainage  area  -30,3  sq  mi0 

;e,— Wonrecording  gage  July  22,  1948,  to  Sept.  30,  1953 J  recording  gage  thereafter, 
Datum  of  gage  is  661,47  ft  above  mean  sea  level,  datum  of  1929 o 

Stage-discharge  relation, —Defined  by  current-meter  measurements  below  330  cfs. 

Flood  3tage,=-10  ft 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


r 

Date 

Gage 

Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

5 

Feb.  15,  1949 

6.15 

254 

1955 

Jan-.  6,  1955 

4-84 

126 

0 

Apr.  4,  1950 

3  3 

422 

1956 

Apr,  29,  1956 

6065 

225 

i 

July  9,  1951 

7.2 

235 

1957 

Apr.  28,  1957 

7.90 

290 

2 

Apr,  23,  1952 

7.3 

263 

1958 

June  10,  1958 

308 

J 

Mar,  15,  1953 

5.75 

185 

1959 

Feb,  10,  1959 

343 

I 

Mar.  25,  1954 

4. 59 

130 

55. 
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(S-10)  Bice  ditch  near  South  Marion,  Ind 


Location, ~~Lat  40°52' ,  long  87°068,  on  line  between  sees,,  15  and  22,  T=  28  N-, 

Ro  6W»,  on  left  bank  at  upstream  side  of  bridge  on  State  Highway  16,  2  miles 
upstream  from  Big  Slough  Creek,  3  miles  southeast  of  South  Iterion,  and  5  miles 
southeast  of  Rensselaer- 

Drainage  area , =-22,6  sq  mi., 

Gage, -- Nonre cording  gage  DeCo  31,  1948,  to  Aug,  4,  1955;  recording  gage  thereafter. 
Datum  of  gage  is  653-30  ft  above  mean  sea  level,  datum  of  1929, 

Stage-discharge  relation, —Defined  by  current-meter  measurements, 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


ter 

ar 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

49 

Feb.  15,  1949 

9  09 

410 

1955 

June  11,  1955 

10,16 

353 

50 

July  19,  1950 

10,06 

490 

1956 

Apr,  29,  1956 

10,75 

504 

51 

July  9,  1951 

11  43 

610 

1957 

July  13,  1957 

10,36 

458 

52 

June  lkB   1952 

10,80 

556 

1958 

June  13,  1958 

780 

53 

July  5*  1953 

fto65 

374 

1959 

Feb,  10,  1959 

480 

54 

Jane  22,  1954 

8,12 

339 

57. 
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(S-ll)  Big  Slough  Creek  near  Collegeville,  Ind< 


Location,  — Lat  40°53»,  long  87°09s,  in  SWj  NW£  sec.  7,  T,  28  N,,  R.  6  w .,  on  right 

bank  on  downstream  side  of  bridge  on  State  Highway  53,  l£  miles  south  of  College- 
ville, 2|  miles  upstream  from  mouth,  and  2  3A  miles  downstream  from  Bice  ditch 

Drainage  araa  <™84«1  sq  mi. 

Gage ,—Nonre cording  gage  July  28,  1948,  to  Dec  31,  1951,  and  Oct.  1,  1952,  to  Aug,  U, 
1955;  recording  gane  since  Aug.  5,  1955.  Datum  of  gags  is  637-. 75  ft  above 
mean  sea  level,  datum  of  1929, 

Stage=dlscharge  relatione --Defined  by  current=meter  measurements. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


.er 
tr 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

.3 

March    1913 

13^7 

1954 

June  22,  1954 

8,84 

390 

!7 

1927 

12.5 

1955 

June  11,  1955 

12,4 

1,100 

9 

Feb.  15,  1949 

11.26 

880 

1956 

Apr  29,  1956 

13  0 

1,470 

0 

Apr„  4,  1950 

12,3 

1,180 

1957 

July  14,  1957 

12  =  96 

1,470 

1 

July  10,  1951 

13-22 

1,450 

1958 

June  13,  1958 

2,030 

3 

Apr.  1,  1953 

10.50 

655 

1959 

Apr,  28,  1959 

1,030 

59. 
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(S-12)  Carpenter  Creek  at  Egypt,  Ind< 


Location , — Lat  40°52',  long  87°12;,  en  line  between  Sw£  sec.  15  and  NW£  sec,  22, 
T.  28  No,  R„  7  Wo,  on  left  bank  on  downstream  side  of  bridge  on  State  Highway 
l6a  2  3/4  miles  upstream  from  mouth,  and  4  miles  southwest  of  Collegeville 

Drainage  Area  —48  1  sq  mi. 

Gage,— Nonre cording  gage  July  26,  1948,  to  Dec  31,  1951,  and  Oct*  1,  1952,  to  Sept  5; 
1955;  recording  gage  since  Sept,  6,  1955c  Datum  of  gage  is  640„37  above  mean 
sea  level,  datum  of  1929 » 

Stags-discharge  relation.— Defined  by  current-meter  measurements. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


e 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

►9 

Feb  15 £  1949 

10  14 

1,160 

1955 

June  8,  1955 

9  80 

984 

.0 

Apr,  4,  1950 

10,3 

1,300 

1956 

Apr„  29,  1956 

9. 93 

1,040 

>i 

July  9,  1951 

10-92 

1,790 

1957 

July  13,  1957 

9*42 

810 

0 

July  6,  1953 

9,21 

735 

1958 

June  10,  1958 

3,720 

14 

June  22,  1954 

8,95 

655 

1959 

Feb,  10,  1959 

2,690 

61. 
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(1)  Tippecanoe  River  at  Oswego,  Ind , 


Location,—  Lat  41°19J14",  long  B5°k7*SL»,   in  NE$  NEj  secc  14,  T„  33  N.,  R  6  E.,  on 
left  bank  10  ft  downstream  from  dam  at  Tippecanoe  Lake  outlet  in  Oswego,  3 
miles  east  of  Leesburgo 

Drainage  area.— 115  sq  mi. 

Gaga ,  -— Nonrecording  gage  Oct,  1,  1949.,  to  Aug.  11,  1953;  recording  gage  thereafter, 
Datum  of  gage  is  830,00  ft  above  mean  sea  level,  datum  of  1929c 

Stage -discharge  relation, —Defined  by  current-meter  measurements  below  680  cfs  and 
extended  to  1,050  cfs  by  logarithmic  plotting „ 

Remarks. —Peak  discharges  affected  by  natural  storage  in  numerous  lakes  upstream. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


sr 

r 

Da*e 

Gage 
Height 

Discharge 
cfs 

Water 

Year 

Gage 
Date       Height 

Discharge 
cfs 

3 

Kay 

21? 

1943 

9  4 

1,050 

1955 

Oct,  17,  1954   8,65 

700 

0 

Apr, 

8— 10, 

1950 

8.62 

680 

1956 

May  5,  6,  1956   8.08 

450 

1 

Feb 

27,28, 

1951 

430 

1957 

Apr,17-22,1957   7,59 

315 

2 

Jan-, 

30, 

1952 

370 

1958 

Sept.  19— 23,1957 

383 

3 

Mar, 

22,23, 

1953 

6.76 

179 

1959 

Feb.  18,  1959 

548 

4 

Apr  -, 

29,30, 

1953 

7.60 

322 

63. 
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(2)  Kississinewa  River  near  Ridgeville,  Ind 


Location,  — Lat  kCPlV ,  long  8500',  in  SW  1/4  sec  8,  T.  19  N„,  Rr.  14  B.f  on  right 
bank  10  ft  downstream  from  highway  bridge,  0o8  mile  downstream  from  Mud  Creek, 
and  2  miles  east  of  Ridgeville„ 

Drainage  area ,=-130  sq  mi, 

Gage— Nonrecording  gage  Aug,  30,  1946,  to  Oct,  3,  1950;  recording  gage  thereafter. 
Datum  of  gage  is  965*23  ft  above  mean  sea  level,  datum  of  1929, 

Stage -discharge  relation c~<-De£ined  by  current -meter  measurements  below  3,400  cfs- 

Flood  stage.— 10  ft. 

Historical  data.=-Local  residents  stated  that  the  1913  flood  was  secondary  to  a 

flood  in  the  early  1930' s  when  the  river  reached  an  estimated  Btage  of  15*0  ft 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


— 

.er 

X 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

7 

Jan, 

30, 

1947 

1116 

2,490 

1954 

Mar,  30,  1954 

7,80 

1,020 

8 

Jan„ 

1. 

1948 

12.2 

3,480 

1955 

Jan.  6,  1955 

11.16 

2,490 

9 

Jan. 

5, 

1949 

13  ol 

4,560 

1956 

Nov.  16,  1955 

12-79 

4,200 

0 

Feb, 

14, 

1950 

13,4 

4,920 

1957 

June  28,  1957 

Ho  57 

8,830 

1 

Feb, 

21, 

1951 

12.75 

4,200 

1958 

June  10,  1958 

13,900 

2 

Jan, 

26, 

1952 

11=99 

3,250 

1959 

Jan0  21,  1959 

^iao 

3 

Mar, 

4, 

1953 

12,00 

3,250 

65. 
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(3)  Wildcat  Creek  at  Greentown,  IncL 


Location.— Lat  40°27",  long  85°57",  on  line  between  sees,  9  and  10,  T.  23  N„,  R  5  E., 
on  left  bank  at  downstream  side  of  bridge  on  State  Highway  213 >  1»5  miles  south 
of  Greentown, 

Drainage  area „— 162  sq  mi;  172  sq  mi  prior  to  June  5,  1954* 

Gage.— Nonrecording  gage  Feb.  20,  1945,  to  June  4,  1954;  recording  gage  thereafter. 
Prior  to  June  5,  1954*  at  site  2.0  miles  downstream  at  datum  5»34  ft  lower  than 
present  datum.  Datum  of  present  gage  is  809.33  ft  above  mean  sea  level,  datum 
of  1929c 

Stage-discharge  relation.  =-=Defined  by  current-meter  measurements. 

Flood  stage. -=-11  ft  at  both  sites. 

Historical  data.==The  following  statements  appear  in  old  newspapers  for  Kokomo, 
about  9  miles  downstream.  September  1866:   "Wildcat  Creek  raging;  no  train 
for  3  days."  1904:  "Greatest  flood  in  Kokono  history." 

Flood  of  August  1871  reached  a  stage  3  inches  below  that  of  the  1913 
flood  at  a  bridge  l£  miles  downstream  from  present  site  according  to  informa- 
tion by  local  resident  on  the  basis  of  remembered  high-water  marks  made  on  the 
same  tree., 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


er 
r 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

3 

May, 

1943 

15.0 

5,960 

1952 

Mar. 

11, 

1952 

11.52 

2,580 

5 

Apr, 

1,  1945 

10,04 

1,680 

1953 

Mar, 

4, 

1953 

10.34 

1,810 

6 

Oct, 

2,  1945 

9-94 

1,640 

1954 

Apr , 

12, 

1954 

5  63 

450 

7 

Apr, 

30,  1947 

10,94 

2,140 

1955 

Jan, 

7, 

1955 

10.00 

1,650 

8 

Mar 

22,  1948 

11.63 

2,670 

1956 

May 

28, 

1956 

9,97 

1,650 

9 

Jan, 

19,  1949 

13  19 

4,110 

1957 

June 

29, 

1957 

12.37 

2,260 

0 

Jan. 

4,  1950 

15.3 

6,320 

1958 

June 

10, 

1958 

4,900 

1 

Feb. 

21,  1951 

10  7 

2,020 

1959 

Feb.. 

10, 

1959 

5,390 
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(4)  Cicero  Creek  near  Arcadia 


Location,™ Lat  40°11» ,  long  86°00! ,  on  line  between  sees.  18  and  19,  T„  20  M  , 
R.  5E.,  on  left  bank  on  downstream  3ide  of  county  bridge,.  l£  miles  east  of 
Arcadia,  Hamilton  County,  and  5  miles  upstream  from  little  Cicero  Creek, 

Drainage  area „— 131  sq0  mi, 

Gage,—  Hater-stage  recorder.  Datum  of  gage  is  815 .  12  ft  above  mean  sea  level 

Datum  of  1929=  Prior  to  Dec.  7,  1955,  wire-weight  gage  at  same  site  and  datum. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 
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sar 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

>55 

July  16,  1955 

1,280 

1958 

June  15,  1958 

2,740 

)5i 

July  21,  1956 

1,540 

1959 

Feb,  11,  1959 

2,170 

>57 

June  29,  1957 

6,720 
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(5)  Fall  Creek  near  Fortville, -Ind, 


Location, --Lat  39°57'>15",  long  85°52»  05",  in  sec  5,  T,  17  No,  R.  6  E .  ,  on  right 
bank  at  downstream  side  of  bridge  on  State  Highway  238*  1  mile  dov.nstream  from 
Lick  Creek  and  2  miles  northwest  of  Fortville 0 

Drainage  area 0=~ 172  sq  mio 

Gage»==Nonrecording  gage  July  1,  1941,  to  June  26,  1942;  record  ijig  gage  thereafter 
Datum  of  gage~is  787 -,43  ft  above  mean  sea  level,  datum  of  1929  (levels  by 
Indianapolis  "Water  Co.,) 

Stage-discharge  relation, —Defined  by  current-meter  measurements, 
Flood  stage  .,=—6  fto 

Historical  data  .—Flood  of  1913  reached  a  stage  op  about  12  feet  according  to  informa- 
tion from  local  resident. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


sar 

•   Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

?42 

Mar,  17, 

1942 

7-39 

2,460 

1951 

Feb, 

22, 

1951 

8„36 

4,250 

943 

May     18  9 

1943 

9,77 

8,240 

1952 

Jan, 

27, 

1952 

6,85 

2.000 

944 

Apr.,   12, 

1944 

80  79 

5,000 

1953 

July 

6, 

1953 

8,14 

3,850 

945 

June  17, 

1945 

6.86 

1,840 

1954 

Mar, 

30, 

1954 

5o50 

1,140 

946 

Oct,      29 

1945 

6,87 

1,840 

1955 

Jan, 

6, 

1955 

5  76 

1,260 

947 

July  15, 

1947 

7»25 

2,250 

1956 

Feb., 

28, 

1956 

7,93 

3,130 

948 

Mar a  24, 

1948 

7,97 

3,500 

1957 

June 

29, 

1957 

8,12 

3,430 

949 

Jan,  19, 

1949 

7d9 

2,350 

1958 

June 

14, 

1958 

5,040 

950 

Jan ,     U, 

1950 

8,85 

5,150 

1959 

Feb, 

10, 

1959 

3,010 

71. 
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(6)  Eagl:     c  at  Indianapolis,,  I 


Location,— Lat  39  46*40%  long  86°15  02!',  in  NW£  sec,  6,  T,  15  N ,  „  B.  3  E»,  on  right 
bank  at  downstream  side  of  bridge  on  Lynnhurst  Drive *  3,0  miles  upstream  from 
Little  Eagle  Creek,  5-0  miles  west  of  Konument  Circle  in  Indianapolis,  and 
6c 7  miles  upstream  from  mouth „ 

Drainage  area— 179  sq  mi. 

Gage,— Nonrecording  gage  Nov.  18,  1933,  to  Mar,  19,  1939;  recording  gage  thereafter 
Datum  of  gage  is  706^21  ft  above  mean  sea  level,  datum  of  1929. 

Stage-discharge  relation , —Defined  by  current-meter  measurements  below  9,000  cfs  and 
extended  above  on  basis  of  a  combined  currents-meter  measurement  and  slop-area 
measurement.:.  High -water  relation  has  shown  a  tendency  to  shift.  Discharge 
shown  for  the  1913  flood  is  an  approximate  value  based  on  logarithmic  extension 
of  an  early  rating  curve  above  9S000  cfs. 

Historical  data.  —The  following  information  was  obtained  from  a  report  on  "Eagle 

Creek  at  Indianapolis,  Channel  Improvements  for  Flood  Control"  by  Indiana  Flood 
Control  and  Water  Resources  Commission,  dated  February  1956:  "Investigations 
on  past  flooding  by  searching  newspaper  files,  interrogation  of  local  persons, 
examination  of  old  surveys  and  reports,  indicate  that  flooding  occurred  along 
Eagle  Creek  in  1875,  1904,  1913,  1915,  1926,  1933 ,  1943,  and  1950,  Newspaper 
accounts  of  flooding  on  other  streams  in  the  Indianapolis  area  indicate  that 
flooding  probably  also  occurred  in  1847,  1858,  1866,  1882  and  1883," 

"It  is  probable  that  the  floods  of  1875  and  1904  were  among  the  greater 
floods  on  the  stream.  Newspaper  accounts  and  weather  records  indicate  that  the 
flood  of  July  1915  was  nearly  as  great  a3  that  of  March  1913." 

Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


rater 
ear 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

913 

March   1913 

16,  C 

11 ' "■  ■  ' 

19,000 

1947 

Jan. 

30, 

1947 

8,47 

3,370 

938 

April   1938 

14  5 

" 

1948 

Apr, 

65 

1948 

12,26 

9,550 

939 

Mar,.  12,  1939 

10- .6 

6,610 

1949 

Jan,. 

19, 

1949 

11.86 

7,250 

940 

Mar,  3*  1940 

6,30 

1,860 

1950 

Jan 

4? 

1950 

13  03 

8,670 

941 

June  12, j  1941 

5o77 

1,470 

1951 

Feb, 

21> 

1951 

8  47 

3,950 

942 

Feb-  7,  1942 

8,66 

4,120 

1952 

Jan. 

27, 

1952 

9,38 

5,520 

943 

May  11,  1943 

1217 

9,660 

1953 

Mar 

5, 

1953 

9-34 

4S920 

944 

Apr,  11,  1944 

10  61 

6,610 

1954 

Apr 

6, 

1954 

6,41 

2,250 

945 

Mar  31,  1945 

8  80 

4,230 

1955 

July 

16, 

1955 

6,91 

2.650 

946 

May  17,  1946 

8.99 

3,860 

1956 
1957 
1958 

1959 

May 
June 
Aug- 
Jan- 

28, 

28. 

8, 

22, 

1956 
1957 
1958 
1959 

13.62 
16, 38 

9,920 

28.800 

8,560 

6,290 

73 , 
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(7)  Youngs  Creek  near  Edinburgh  Ind, 


Locatlono— Lat  39  25*08",  long  86  0018%  in  SW  1/4  sec.  5,  Ts  11  MOJ  R..  5  e.„  on 
left  bank,  on  upstream  side  of  highway  bridge  half  a  mile  southwest  of  Amity, 
2  miles  upstream  from  mouth,  and  5  miles  northwest  of  Edinburg, 

Drainage  area  —109  sq  mi. 

Gage ,~-Nonre cording  gage  Dec*  7?  1942,  to  June  29,  1955;  recording  gage  thereafter 
Datum  of  gage  is  670,20  ft  above  mean  sea  level,  datum  of  1929. 

Stage^discharge  relation .—Defined  by  current -meter  measurements  below  7,000  cfs 
and  by  contracted«opening  measurement  at  10^700  cfs. 

Flood  stage*— 7  ft. 


Feak  Stages  and  Instantaneous  Annual  Peak  Discharge 


ter 
ar 

43 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

Mar 

19,  1943 

10,40 

3*700 

1952 

Janc  27,  1952 

13  4 

10,700 

44 

Apr, 

11,  1944 

11,00 

4,290 

1953 

Mar,  4p  1953 

8,37 

2,080 

45 

Mar, 

6,  1945 

11-00 

4,290 

1954 

Jan,  27,  1954 

3  27 

443 

46 

May 

16,  1946 

9  0 

2,510 

1955 

May  28 ,  1955 

6-2 

1,110 

W 

June 

2,  1947 

11,12 

4,390 

1956 

Nov,  16,  1955 

12.20 

7,790 

*8 

Mar- 

27,  1948 

7.68 

1,650 

1957 

July  5,  1957 

11.62 

6,510 

fl 

Jan 

5.  1949 

11  9 

5,190 

1958 

June  11,  1958 

4,350 

30 

Jan 

4,  1950 

10-8 

4,090 

1959 

Jan,  21,  1959 

6,270 

51 

Jan., 

15,  1951 

9.15 

30 

75< 
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(8)  Blue  River  at  Carthage,  Ind„ 


L0Cat^  "Lat  ?9°^>  W  35°34»,  in  sec.  18,  T  15  N„  R  9  E.,  on  right  bank 
500  ft  upstream  from  highway  bridge,  half  a  mile  west  of  Carthage,  and  2  1 
miles  downstream  from  Three  Mile  Creeko 

Drainage  area  „  —187  sq  rai 

Gage  .^Nonrecording  gage  Of  -  Hs  1950,  to  July  18,  1951  j  recording  gage  thereafter 
Prior  to  July  19,  1951?  at  bridge  500  ft  downstream.  Datum  of  gape  is  859  33  ft 
above  mean  sea  level s  datum  of  1929, 

Stage-discharge  relation „ —Defined  by  current-meter  measurements, 
Flood  stage  -=7  ft. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


iter 
;ar 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

>49 

Jan.      5,  1949 

10,6 

5*750 

1955 

Jan,     6,  1955 

597 

1,290 

>a 

Feb.    21,  1951 

11,2 

6,650 

1956 

Nov.  16,  1955 

11,52 

5,800 

>52 

Jan,  27,  1952 

11  02 

6,350 

1957 

June  18,  1957 

977 

3,900 

•53 

Mar,     4,  1953 

9,17 

38580 

1958 

June  14*  1958 

7,020 

% 

Apr,     5,  1954 

10  02 

4,850 

1959 

Jan,  21,  1959 

8,340 

J 
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(9)  Sand  Creek  near  Brewersville,  Ind, 


Location.— Lat  39°05805",  long  85°39v30»,  in  NW$  sec  5,  To  7  N.  ,  R.  8  E,,  on  left 

bank  at  downstream  side  of  county  highway  bridge,  2\  mileE  west  of  Brewersville, 
and  5  2  miles  upstream  from  Bear  Creeko 

Drainage  area-— 156  sq  mi;  153  sq  mi  prior  to  Oct,  6,  1952„ 

Gage.==Nonrecording  gage  Febo  11,  1948,  to  Oct.  5,  1952,  at  bridge  1„7  miles  upstream 
at  datum  approximately  8  ft  higher.  Recording  gage  since  Oct,  6,  1952,  at 
present  site.  Altitude  of  present  gage  is  630  ft  (by  altimeter). 

Stage-discharge  relation.,— Defined  by  current-meter  measurements  at  former  site  and 
by  gage-height  relationship  with  former  site  at  present  location 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


ter 
ar 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

IS 

Mar     27;  1948 

17-5 

9,980 

1954 

June     1,  1954 

5,75 

1,240 

49 

Jan,     5 ,  1949 

19  0 

12,3.00 

1955 

Feb     27,  1955 

11  42 

u  300 

50 

Jan      4 *  1950 

19  2 

12,400 

1956 

May  28,  1956 

15  45 

7,560 

51 

Nov,  20,  1950 

18c  4 

11,100 

1957 

Apr,     4,  1957 

16.33 

8,480 

52 

Jan     26,  1952 

13  4 

5,780 

1958 

July  22,  1958 

7,150 

53 

Mar      4,  1953 

10.19 

3,460 

1959 

Jan.   21,  1959 

19,900 
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(10)  North  Fork  Salt  Creek  near  Belmont,  Ind- 


Location,  —  Lat  39  09600»8  long  86  20'14",  in  NW  l/4  aeco  5,  T,  8  No,  i     2  B«,  on 

right  bank  15  ft  downstream  from  bridge  on  State  Highway  46  ,  100  ft  upstream 
from  Schooner  Creek,  0,7  mile  northeast  of  Belmont,  6  1/2  miles  Upstream  from 
Brummett Creek ^  and  20  miles  upstream  from  mouth, 

Drainage  area,—120  sq  mis  includes  that  of  Schooner  Creek, 

Gage  .==Nonrecording  gage  Apr,  4,  1946,  to  Oct,  8,  1951;  recording  gage  thereafter 
Altitude  of  gape  is  546  ft  (from  topographic  map) 

Stage-discharge  relation „ —Defined  by  current-meter  measurements  below  9,800  cfs. 

Discharge  adjusted  for  rate  of  change  of  stage  above  7  ft:  Only  annual  maximums 
adjusted  prior  to  installation  of  recording  gage, 

Flood  stage.,— 16  ft 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


;er 

IT 

Date 

Gage 
Height 

Discharge 
cfs 

V/ater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

3 

March 

1913 

25,7 

- 

1953 

Mar   4, 

1953 

17o76 

2;780 

,6 

May  16, 

1946 

20,3 

5,910 

1954 

Jan,  27, 

1954 

9. 38 

825 

7 

June  2, 

1947 

212 

10*100 

1955 

Mar„  21, 

1955 

15-91 

2  220 

3 

Mar,  27, 

1948 

18  o0 

3,010 

1956 

May  28, 

1956 

18.12 

3,030 

9 

Jan,,  5j> 

1949 

22  2 

13,300 

1957 

Apr.  4, 

1957 

19  c  92 

6,340 

so 

Jan ,  4, 

1950 

21  7 

11,600 

1958 

June  14, 

1958 

5,920 

ii 

Feb,  21, 

1951 

19-53 

5*100 

1959 

Jan.  21, 

1959 

000 

■  2 

May  24, 

1952 

22  55 

15,200 

31. 
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(11)  Patoka  River  at  Jasper ,  Ind . 


Location.— Lat  38024'49",  long  86c52<36%  in  SE£  sec,  20,  T.  1  S.  ,  R,  4  W,,  on  left 
bank0  0  3  mile  upstream  from  unnamed  outlet  of  Jasper  Lake?  1  0  mile  downstream 
from  Coon  Seitz  bridge P  1.2  miles  downstream  from  Beaver  Creek^  and  3,3  miles 
northeast  of  Jasper, 

Drainage  are?.  .—257  sq  mi;  270  3q  mi  at  former  site 

Gage  —  Nonrecording  gage  Nov.  20  s  1947,  to  Sept,  17,  1956;  recording  gage  thereafter 
Prior  to  Sept  18,  1956,  at  site  5»6  miles  downstream  at  datum  0,34  ft  lower; 
datum  of  present  gape  is  446,19  ft  above  mean  sea  level,  datum  of  1929, 

Stage  discharge  relation,,  —  Defined  by  current-meter  measurements  below  5,000  cfs  at 
former  site  and  below  1,100  cfs  for  present  site0 

Flood  stage  --14  ft;  9  ft  at  former  site. 

Historical  data, —Flood  of  March  1913  is  maximum  stage  known.  Maximum  stage  at 

present  site  for  period  1925-57,  20  ft  in  1925  (information  from  local  resident). 

Remarks „— ■  Flow  slightly  regulated  by  Beaver  Creek  Reservoir,  whose  outlet  enters  the 
Patoka  River  1*2  miles  upstream  from  the  gage;  peak  discharges  not  materially 
affected. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


sr 

r 

Date 

Gage 
Height 

Discharge 
cfs 

Water 

Year 

Date 

Gage 
Height 

Discharge 
cfs 

3 

March 

1913 

15-9 

16,000 

1953 

Mar, 

8,  1953 

7  90 

1,640 

7 

January 

1937 

14,8 

12,100 

1954 

Mar, 

2,  1954 

- 

950 

S 

Apr,  15, 

1948 

11-57 

4,920 

1955 

Mar. 

3,  1955 

9.80 

2,940 

9 

Jan*  28, 

1949 

11,13 

4,220 

1956 

Feb, 

29,  1956 

9»98 

3,100 

0 

Jan,  7, 

1950 

12.37 

6,300 

1957 

May 

25,  1957 

17,87 

6,900 

1 

Mar,  21, 

1951 

11.46 

4,760 

1958 

Dec. 

22,  1957 

4,250 

2 

Mar.  14, 

1952 

10.78 

3,880 

1959 

Jan-. 

24,  1959 

9,150 
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(12)  Busseron  Greek  near  Carlisle,  Ind, 


Location.—  Lat  38°58  30%  long  87°25tl35"J  in  w£  sec-,  17,  T,  6  JC,  R.  9  W-.,  on  right 

bank  10  ft  downstream  from  bridge  on  State  Highway  58,  l|  miles  northwest  of  Car- 
lisle, and  6  3/4  miles  upstream  from  mouth. 

Drainage  area  =  —228  sq  raia 

Gage.,—Nonrecording  gage  Oct.  15 ,  1943*  to  Nov,  7,  1950 j  recording  gage  thereafter 
Datum  of  gage  is  425^36  ft  above  mean  sea  level  (State  Highway  Department  of 
Indiana  bench  mark). 

Stage-discharge  relation „ -'-Defined  by  current=meter  measurements  below  4,500  cfs  and 
extended  above  by  logarithmic  plotting. 

Flood  stage. —12  ft. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


iter 

ear 

Date 

Gage 

Height 

Discharge 
cfs 

Kater 
Year 

Date 

Gage 
Height 

Discharge 

cfs 

?44 

Apr,  12, 

1944 

16  96 

4,700 

1952 

Kar,  11, 

1952 

16.17 

4,070 

?45 

Apr,     2, 

1945 

17,60 

5,500 

1953 

Mar,     4, 

1953 

16,04 

3,890 

)46 

May    20, 

1946 

14  90 

2,900 

1954 

Aug,     4, 

1954 

6  31 

430 

>47 

June     2, 

1947 

14,60 

2,720 

1955 

Apr.  13, 

1955 

13.13 

2,040 

)ue 

Jan.     3, 

1948 

15-15 

3A00 

1956 

June  22, 

1956 

16.12 

3,980 

>49 

Jan,  20, 

1949 

16.3 

4,200 

1957 

May     23, 

1957 

17,61 

5,200 

?50 

Jan.,     5, 

1950 

20.05 

8,800 

1958 

Dec     21 p 

1957 

5,400 

?51 

Feb,   21, 

1951 

14,75 

2,900 

1959 

Jan,   22, 

1959 

3,100 

85. 
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(13)  East  Fork  V/hitewater  River  at  Richmond,  Ind; 


Location. —Lat  39°4Sl24",  long  84°54«26",  in  3E  1/4  sec,  7,  T.  13  R.,  R.  1  W.,  on 
left  bank  50  ft  downstream  from  highway  bridge,  three-quarters  of  a  mile  south 
of  Richmond,  and  2  miles  upstream  from  Short  Creek. 

Drainage  area. — 123  sq  mi, 

Ga^s . --Nonrecording  gags  Apr*  27,  1949,  to  July  26,  1949;  recording  gage  thereafter. 

Datum  of  gage  is  354  01  ft  above  mean  sea  level,  datum  of  1929  (levels  by  Indiana 
Flood  Control  and  Water  Resources  Commission). 

Stase-discharge  relation. —Defined  by  current-meter  measurements  below  5, ICO  cfs  and 
by  slope-area  measurement  at  13,500  cfs. 

Flood  stage, =-10  ft. 

'Historical  data.— Flood  of  September  1866  was  reported  by  the  Indianapolis  Journal  to 
be  higher  than  ever  before  known.  Flood  of  March  1913  ia  the  maximum  stage  known 
according  to  information  by  local  residents. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


J5>tt?  V* 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

.913 

March    1913 

15.0 

1955 

Feb  21,  1955 

6  07 

2,540 

■950 

Jan  15,  1950 

12,49 

13s 500 

1956 

Nov  16,  1955 

10,70 

8,200 

951 

Nov,  20,  1950 

10.82 

9,250 

1957 

June  28,  1957 

10-54 

7,800 

■952 

Jan,  26,  1952 

10,66 

9,250 

1958 

Aug,  2,  1953 

10,400 

-953 

May  22,  1953 

6.53 

3,160 

1959 

Jan.  21,  1959 

14,100 

-95i 

Mar,  30,  1954 

3.86 

1,160 
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(14)  Silver  Greek  near  Sellersburg 

tion.—Lat  38°22i15",  long  85°43'35",  in  SW£  Lot  68,  Clark  Military  Grant, 
upstream  side  of  Straws  Kill  3rldge  oft  V,atson  Road,  0.3  mile  downstream  from 
Pleasant  Hun,  2,4  miles  southeast  of  Sellersburg,  and  11, S  miles  upstream  from 
mouth „ 

iage  area o— 188  sq„  mi. 

:.-»Wire -weight  gage  read  vwice  daily. 


Peak  Stapes  and  Instantaneous  Annual  Peak  Discharge 


it 


Date 


19; 


Feb,  22,  1955 

Feb,  2,  1956 

23,  1957 


Gage 
height 


Dis 


ks  220 
5,250 
6,250 


Water 

Year 


Date 


Gage    Discharge 
Height      cfs 


1953 
1959 


Nov,  19,  1957 
Jan  22,  1959 


19.600 
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(15)  Big  Indian  Creek  near  Corydon,  Ind. 

Location.  -=Lat  38°l6»35",  long  86o06*35",  in  S3  1/4  sec,  6,  T.  3  S.,  R,  4  E.,  on 
upstream  side  of  bridge  on  State  Highway  335,  0,6  mile  upstream  from  iiaccoon 
Branch  and  4  l/2  miles  north  of  Corydon, 

Drainage  area,~«129  sq  mi. 

Gage     irecording  gage  Oct.  12,  1943,  to  Dec,  89  1948;  recording  gage  thereafter, 
Datum  of  gage  is  577.12  ft  above  mean  sea  level,  datum  of  1929, 

Stage=ti:.seharge  relation,'— Defined  by  current -meter  measurements  below  6,600  cfs  and 
tended  above  by  logarithmic  plotting. 

Historical  data,-— Flood  of  Mar,  19,  1943,  is  the  maximum  known  at  Corydon  since 
beginning  of  knowledge  in  1815  -. 
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Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

19, 

22 

17,000 

1952 

Mar,  11 , 

1952 

18.76 

10,400 

Apr  . 

11, 

194** 

14  0 

4;  860 

1953 

Mar,  3, 

1953 

14,57 

5.400 

Mar, 

6, 

19,2 

000 

1954 

'Sept<,20, 

1954 

14  04 

4,800 

14, 

15o0 

5:910 

1955 

Mar,  15, 

1955 

12,98 

3,900 

2, 

1947 

14,3 

5,160 

1956 

Feb,  25, 

1956 

16>25 

7,180 

Apr, 

12, 

19,3 

11,100 

1957 

Feb,  10, 

1957 

16.32 

7,300 

15, 

16,20 

7,180 

1958 

Nov.  19, 

1957 

6,590 

Feb, 

9, 

1950 

16  ,.77 

7,900 

1959 

Jan.  21, 

1959 

23,800 

Jan. 

14, 

1951 

15.9 

6 ,  320 

91. 
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(16)  Kankakee  River  near  North  Liberty,  I 


Location,  —Lat  4I°33C50",  long  86c29°50",  on  line  between  sees.  11  and  .23,  To  36  K-., 
5n  left  bank  at  downstream  side  of  bridge  on  5t,  Joseph  County  highway 
cad."  4  miles  northwest  of  Forth  Liberty. 

Drainzge  area      3q  mi. 

Nonrecording  gage  Jan.  12,  1951,  to  June  25,  1956;  recording  gage  thereafter, 
Datum  of  gage  is  680.04  ft  at ova  nean  sea  level,  datum  of  1929  (levels  by  Indiana 
Flood  Control  and  F'ater  Resources  Commission). 

Stage-disc      relation. -=Relat ion  affected  by  varying  amounts  of  backwater  caused 

by  return  flow  from  overbank  storage.   Freouent  current-meter  measurements  neces- 
ry  to  define  relationship  during  this  period. 
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ater 
ear 

Date 

Gage 

Height 

rge 

cfs 

er 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

.951 

. 

6.15 

520 

3-956 

Apr,  30,  1956 

6,92 

660 

952 

rvo   14,  1951 

6.97 

680 

1957 

r.  27,  1957 

6.90 

660 

953 

ar.  16,  1953 

4.42 

260 

1958 

Dec.  20,  1957 

560 

954 

tr    26,  1954 

800 

1959 

Kar.  27,  1959 

560 

955 

Oct,  10,  1954 

8.64 

686 
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(17)  Singleton  ditch  at  Schneider,  Ind. 

Location,,-- Lat  41°12!!44",  long  87°26W,  on  line  between  IIS  1/4  sec.  21  and  NW 

1/4  sec,  22,  T„  32  N.,  U„  9  W.,  on  left  bank  15  ft  upstream  from  bridge  on  U.  S, 
Highway  41,  half  a  mile  upstream  from  Bruce  ditch,  1  i/2  miles  downstream  from 
Cedar  Creek,  and  1  2/3  miles  north  of  Schneider., 

Drainage  area. =-122  sq  mi. 

Gage,=-konrecording  gage  July  23,  1948,  to  Aug.  13,  1951;  recording  page  thereafter. 
Prior  to  Oct.  1,  1949,  at  datum  2.00  ft  higher.  Datum  of  present  gage  is  623  67 
ft  above  mean  sea  level,  datum  of  1929. 

Stage-discharge  relation.— Defined  by  current-meter  measurements,  Dredging  in  1950 
""  and  subsequent  floods  and  channel  deterioration  have  materially  affected  the 
stage-discharge  relation. 
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Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1949 

Feb-  15,  1949 

- 

550 

1955 

Oct. 

11,  1954 

10.10 

953 

1950 

Apro  10,  1950 

- 

1,100 

1956 

Feb. 

25,  1956 

9-62 

888 

1951 

Feb,  19,  1951 

8,50 

841 

1957 

Apr, 

28,  1957 

10.27 

979 

1952 

June  15*  1952 

9.82 

1,010 

1958 

714 

1953 

8.39  ' 

812 

1959 

Feb. 

14,  1959 

992 

1954 

Mar„  25,  1954 

9o04 

810 
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(18)  Deep  River  at  Lake  George  Outlet  at  Hobart,  Ind. 

Location*— Lat  41°32»10»,  long  ST^-S^",  in  NW  1/4  sec,  32,  T,  36  N.,  R0  7  W>,  on 
left  bank  at  upstream  side  of  highway  bridge,  300  ft  upstream  from  Duck  Creekt 
and  400  ft  downstream  from  Lake  George  Dam. 

Drainage  area, -=125  sq  mi. 

Gage o— Nonre cording  gage  Apr.  4,  1947,  to  July  29,  1952;  recording  gage  thereafter. 
Prior  to  July  21,  1955,  at  site  400  ft  upstream  at  datum  11.80  ft  higher  than 
present  datum.  Datum  of  present  gage  is  588=17  ft  above  mean  sea  level,  datum 
of  1929  (levels  by  Indiana  Flood  Control  and  Iv'ater  Resources  Commission). 

Stage=discharge  relation. —Defined  by  current-meter  measurements  below  3,300  cfs. 


Peak  Stages  and  Instantaneous  Annual  Peak  Discharge 


feiter 

'ear 

Date 

Gage 

Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

947 

Apr.  6,  1947 

5.41 

2,410 

1954 

Mar-  26,  1954 

4.55 

1,440 

948 

May  11,  1948 

5,86 

2,740 

1955 

Oct„  11,  1954 

7.68 

3,880 

949 

Feb,  14,  1949 

3o50 

620 

1956 

Kay  11,  1956 

11.15 

1,320 

950 

Dec.  22,  1949 

5 -.35 

2,390 

1957 

July  14,  1957 

12o35 

1,650 

951 

May  11,  1951 

4,52 

1,440 

1958 

June  10,  1958 

720 

952 

Nov.  14,  1951 

4.41 

1,340 

1959 

July  24,  1959 

1,970 

953 

Mar.  16,  1953 

3.86 

912 

97. 
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(19)  Pigeon  Creek  at  Hogback  Lake  Outlet,  near  Angola*  Ind, 


Location.,— Lat  41037'24!',  long  8:'°05844",  in  KE  l/k  NW  l/4  sec.  36,  T.  37  N.,  R.  12  E., 
on  right  bank  200  ft  north  of  lake  outlet,  2  miles  southeast  of  Flint,  and  5.1 
miles  west  of  Angola, 

Drainage  area. —102  sq  ml*  105  sq  mi  prior  to  October  1947. 

Gage. — Nonre cording  gage  Oct.  16,  1945,  to  Aug.  3,  1953;  recording  gage  thereafter. 

Prior  to  Oct.  1,  1947,  at  site  1  l/2  miles  downstream  at  different  datum.  Oct.  1, 
1947,  to  Aug.  3,  1953,  at  site  600  ft  downstream  at  present  datum.  Datum  of 
present  gage  is  940.00  ft  above  mean  sea  level,  datum  of  1929. 

Stags=<iis  charge  relation. —Defined  by  current-meter  measurements  below  240  cfs  at 
former  site  and  by  current°meter  measurements  at  present  site. 
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Water 
Tear 

Date 

Gage 
Height 

Discharge 
cfs 

Water 
Year 

Date 

Gage 
Height 

Discharge 

Cf3 

1946 

Feb. 

19,  1946 

- 

220 

1953 

Mar.  19,  1953 

9o30 

122 

1947 

Apr. 

24,  1947 

10.71 

458 

1954 

Kar.  30,  1954 

11.31 

317 

194£ 

liar. 

2,  1948 

11.79 

355 

1955 

Oct.  17,  1954 

11  =  54 

339 

1949 

Feb. 

19,  1949 

11.93 

366 

1956 

fey   4,  1956 

13*39 

548 

1950 

Apr. 

8,  1950 

14.95 

744 

1957 

Apr.  14,  1957 

11.29 

317 

1951 

Feb. 

24,  1951 

12.50 

448 

1958 

Sept. 21,  1957 

274 

1952 

Jan, 

21,  1952 

11 085 

370 

1959 

Feb, 17-19, 1959 

442 
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(20)  Laughery  Creek  r.ear  Farmers  ftetreat,  Ind, 


Location.— Lat  38o57<05«,  leng  85°Q4«22«,  in  sec.  2,  To  4  IL,   R„  3  W.,  on  right  bank 
2  miles  southeast  of  Farmers  Retreat  and  3  3/4  miles  downstream  from  Rear  Creek, 

Drainage  area..—  248  sq  mi. 

Gageo=^Ionrecording  gage  0cto  3,  1940,  to  Apr.  15,  1941;   recording  gage  thereafter. 
Altitude  of  gage  is  526  ft  (by  barometer). 

Stage-discharge  relation. —Defined  by  eurrent-^neter  measurements  below  14,000  cfs. 

Flood  stage <,=13  ft. 

Historical  data* —Flood  of  1897  reached  a  stage  of  about  18  feet  and  is  the  highest 
known  flood,  from  information  by  local  residents <, 
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Water 

pfear 

Date 

Gage 
Height 

Discharge 
cfs 

Viiater 
Year 

Date 

Gage 
Height 

Discharge 
cfs 

1941 

June  9 

or 

10,1941 

9^62 

3,960 

1951 

Jan, 

3, 

1951 

13,50 

9,660 

J1942 

Apr, 

9, 

1942 

13o91 

10,800 

1952 

Mar., 

10, 

1952 

13,46 

9,660 

1943 

Mar0 

19, 

1943 

14c  50 

12,900 

1953 

May 

17, 

1953 

9,44 

3,960 

1944 

Apr, 

11, 

1944 

13ol6 

8,880 

1954 

May 

3, 

1954 

3  99 

640 

1945 

Mar, 

6, 

1945 

15, 5A 

17,000 

1955 

Mar, 

21, 

1955 

13-88 

10,800 

1946 

Feb, 

13, 

1946 

12,78 

7,980 

1956 

May 

28  8 

1956 

14,45 

12 t 500 

1947 

May 

25, 

19^7 

14,62 

13,300 

1957 

July 

5, 

1957 

16,15 

20,200 

1948 

Apr, 

12, 

1948 

13.01 

8,410 

1958 

July 

22, 

1958 

17.000 

1949 

Jan, 

24, 

1949 

15o23 

15,700 

1959 

Jan, 

21, 

1959 

47,800 

3, 

1950 

03 

1/  ,800 
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Table  IV 
The  Fredicted  Annual  Instantaneous  Peak  Discharge  ^rom 
Flood  Frequency  Analysis 


Watershed 
No 

Predict 

;ing  Annual  Instantaneous 

Peak  Discharge 

25-yrs 

50-yrs      75=yrs 

100=yrs 

cfs. 

cf»o        cfs. 

cfs. 

S~l 

11800 

13800       U700 

15900 

S~2 

12900 

15000      16000 

17200 

S~3 

20>00 

2LCQ0                 25400 

27000 

S=A 

3300 

3750       3940 

4200 

S=5 

2950 

3420       3600 

3860 

S=6 

3300 

3800       4000 

4300 

S»7 

1630 

1800       1880 

1980 

s~e 

2100 

2350       2450 

26OO 

S=9 

465 

520        540 

570 

S~10 

835 

925        960 

1010 

S-ll 

2440 

2770       2900 

3080 

S-12 

4160 

4800       5100 

5500 

1 

880 

1000       1040 

1110 

2 

14000 

16500      17400 

18800 

3 

7400 

8600       9000 

9700 

4 

6800 

7800       8300 

8900 

5 

7700 

8900       9300 

10000 

6 

17000 

19700       21000 

22300 

7 

10500 

12400       13000 

I4OOO 

8 

10800 

12200       12800 

13600 

9 

21500 

21,000       25200 

26800 

10 

19000 

22000       23300 

25200 

11 

13500 

15600       16500 

17700 

12 

8700 

9900       10300 

11000 

13 

19800 

22800       24200 

26000 

14 

13300 

15100       15900 

17000 

15 

22300 

25700       27300 

29400 

16 

1020 

1130        1170 

1230 

17 

1270 

1370        1410 

1470 

18 

4800 

5500        5700 

6200 

19 

760 

860         900 

960 

20 

28600 

33000      35COO 

37300 
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7c  GEOMORPHOLOGICAL  STUDY 
The  Purpose  of  Geomorphological  Study 

After  a  study  of  the  flood  frequency  analysis  and  obtaining  the 
25  years  flood  for  seal!  gaged  watersheds ,  a  3tudy  was  made  to  estimate 
the  flood  from  ungaged  watersheds „  This  study  investigates  the  geomor- 
phological factors  of  each  watershed  and  attempts  to  correlate  them  to 
the  flood's  flows c  If  a  relationship  can  be  obtained  between  the  geo= 
morphological  factors  describing  the  watershed  and  its  corresponding 
25  year  flood,  then  the  flood  discharge  can  be  estimated  by  simply 
measuring  the  geomorphological  factors  from  a  topographic  map. 

The  geomorph.©  logical  3tudy  is  based  on  the  assumption  that  region 
is  homogeneous.  That  is  to  say  that  the  3torm  pattern  above  and  the 
geological  conditions  below,  are  homogeneous 0  It  is  assumed  that  the 
geomorphological  factors  are  the  saain  factors  that  control  the  rate  and 
amount  of  runoff.  The  first  geomorphological  factor  put  into  considera= 
tion  is  the  drainage  area   It  is  the  projected  area  of  watershed p  also 
called  the  catchment  area.  Obviously s   the  bigger  the  area  of  catchment 
the  larger  the  amount  of  runoff.  However^  the  rats  of  runoff  is  largely 
dependant  on  the  slope  of  the  land,  the  drainage  density,  and  the  slope 
of  the  streams,  A  atesp  land  slope  with  high  drainage  density  and  larger 
slope  of  stream  will  give  a  higher  rate  of  runoff  than  those  which  have 
smaller  -slope  and  lower  drainage  density*  The  shape  of  the  watershed  is 
also  a  factor  af  fee  Ling  Liis  runoff  <,  An  elongated  watershed  will  have 
a  longer  time  of   concentration  than  those  with  round  shapes ,  and  thus 
less  discharge  will  take  place. 

With  a  three-dimensional  concept  in  mind,  the  geomorphological  factors 
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which  can  affect  the  aaount  of  runoff  are  listed  below: 
1)  The  drainage  area 

The  drainage  density 
3)  The  land  slope 

te  main  stream  si 
sd  shs] 

Evaluation  of  Geomorpholo-     actors 

ini'.a—.ii,     ■T.iff.JT  CiJniinm-nMBaa»*»5ManiimiMi«  lYfr  ■■  —  m  mm  inn  i  n,  unuw 

1)  Drainage  Area 

Area  of  the  watersh--    directly  measured  from  topographic  maps 
with  a  planimeter<>  It  i.     :r,sed  in  square  miles. 

2)  Draii  Density  (17) 

Drainag  total  length  of  streams  in  the 

3d  by  Its        a,  that  is  the  length  of  streams  per 
itershedo  It  can  be  expressed  as  flllows: 


Drainage  Density 
D 


(3) 


Where 

L  *  total  length  of  streams 

A  «  drainage  area 


The  length  of  stream  can  be  measured  in  miles  by  map  measurement 

alor  topographic  map-  Since  the  area  is 

expressed  in  square  miles  i;  :ainage  density  can  be  expressed  as 
miles  per  square  miles. 

3)  Land  Slope  (1?) 

Since  the  land  slope  if  changing  from  place  to  place-  in  a  watershed,, 

it  is  hard  to  find  a  quar  e  value  to  represent  the  land  slope  of  a 
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whole  watershed   A  new  parameter  introduced  here  to  replace  the  land 
slope  is  the  mean  relief  of  land,,  The  mean  relief  i3  def ined  as  the  total 
volume  of  land  mass  above  the  gaging  station  of  a  waterskad  divided  by- 
its  projected  area,.  This  can  be  evaluated  quantitatively  by  using  the 
so-called  hypsometric  analysis „  The  hypsometric  curve  developed  by 
Langbain  and  others  (15)  gives  a  dimensionless  relationship  between  the 
horizontal  cross-sectional  drainage  basin  area  and  the  elevation-  The  total 
volume  of  land  mass  can  be  calculated,  and  the  mean  relief  also  determined . 


area-*  /}. 


taosrfb- 


Fig.  U) 

Figure  (4)  shows  the  elevation  and  Its  related  cross-sectional  drainage 
basin  &re&a  The  elevation  can  be  read  from  topographic  mapSj,  and  its 
cro8s«sectional  area  can  be  measured  by  using  the  planimeter,.  By  obtain- 
ing ths  H*,  hg  A,  a,  from  the  topographic  maps  a  hypsometric  curve  can  be 
plotted  in  a  dimanaionleas  form  as  shown  in  Figure  (5). 


Fig.  (5) 


/?*/ai/u»  -4/jffJC 


106. 

From  the  above  hypsometric  curve,  the  total  volume  of  land  mass  and  the 
mean  relief  can  be  determined 0  Since  the  area  under  the  hypsometric  curve 
can  ba  easily  measured,  this  can  be  expressed  as, 


X 


*<&>  •  •<■ 


This  can  be  rewritten  as 


-  < 


The  total  volume  of  land  mass 


,a  Ca*&  -  <*£_/?/ 


The  Esan  relief 


X  = 


H  -     -fl   =  ^J* 


* 


4)  Main  Stream  Slope 

The  slope  of  main  stream  can  be  obtained  from  the  topographic  mapp 
but  it  is  not  the  same  throughout  the  whole  stream,,  Usually *  the  upper 
streams  are  steepsr9  and  the  downstream  reaches  are  flatter.  The  problem 
is  how  to  get  a  mean  slope  to  represent  the  main  stream  slope 0  It  would 
not  be  sufficient  .just  to  find  the  slope  of  a  straight  line  which  connects 
the  uppsr  stream  and  the  downstream  A  method  introduced  by  Taylor  and 
Schwary  (25)  for  determining  raean  slope  of  stream  is  shown  a3  follows: 
(1)  Divide  stream  into  nany  sections,  for  each  increment  of 
stream  length 


Where 


t  ■  travel  time  of  flow  through  the  increment  of  length 

1  ■  increment  of  length 

v,  ■  average  velc-ciiy  of  flow  through  the  increment. 
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(2)  From  Manning" a  formula 


and  for  each  length  of  increment  -~p  R      le  determined 
and  equal  Kj 

thU3, 


V,  =  A",  5 


(3)  total  travel  time 

7~=    fi    +  4    ■/  ?3   +  €4  +    '  '   '    +    *n 

V,  VZ  *3  Vf  ** 

(4)  If  equal  increntents  of  length  are  used 
then, 


5)  Assuming  a  uniform  channel  having  hydraulic  elements 
equivalent  to  the  average  for  the  prototype  channel , 
the  travel  time  would  be 

and  the  velocity 


elope 
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Substitute  this  relation  to 


T-  it  (*&*  *  ^t* 


■^■Sn 


wj 


;■>.  •+■ 


**% 


a 


s 


Y* 


rt 


K« 


TsZ*) 


*«    **    **  *W  *  / 


5)  V.'atorshed  shape  factor 

The  watershed  shape  factor  is  the  factor  that  transforms  the  irregular 
shape  of  a  watershsd  to  a  standard  one  which  is  circular  The  length  of 
longest  waterway  is  used  as  paraiaater  to  determine  the  shape  factor,- 


Fig*  (6) 

Figure  (6)  shows  two  watersheds  with  the  same  area  but  different  shape. 
The  former  is  an  actual  watershed,  and  the  latter  is  an  imaginary  one„ 


.•oo 
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A0  is  measured  along  the  main  stream,  and  A»0«   ia  calculated  from  the  known 
area  of  the  actual  watershed 


area  or  /Qcrg/ii     w<QT£&SH'eo 

*  b'  -    2  /  tP 

by  assuming  that  this  watershed  is  a  circular  one  and  A'O6  is  the  diameter 
of  the  circle,. 

The  shape  factor  is  determined  as      f  «  ...ft  7  — 


Geomorphological  Factors  of  Small  Watoraheda  in  Indiana 

The  following  table  (V)  shows  quantitatively  the  gaomorphological 
factors  of  16  small  watersheds  in  Indiana. 
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Table  V 
Watershed  Characteristics  and  25-Year  Annual  Instantaneous 
Peak  Runoff  cf  16  Watersheds  in  Indiana 


Watershed 

r : , 

25=Year  Annual 
:antaneous 
Peak  Runoff 

Q  Cf3o 

...  Wate_rj 

Jhed    Char 

Acteristics 

Area 
A  sc^tti. 

Mean 

Relief 
H  ft. 

Drainage 
Density 
D  mi/sq.mi. 

Shape  naon  stream 
Factor    Slope, 
f    S  x  lO** 

S=l 

11, SCO 

100 

216 

10.66 

2c  63 

-4 

9.84x  10 

S-2 

12,900 

88,8 

270 

7.38 

300 

20,88 

S=5 

2,950 

78o7 

101 

6,57 

lc75 

9-05 

S=6 

3,300 

62.>9 

110 

8,00 

XqJLa 

21,10 

3-7 

1,630 

93-0 

79 

5,10 

1*47 

8,29 

1 

880 

115 

65,4 

3.35 

lo4l 

2,64 

6 

17,000 

179 

195.2 

7.88 

2,12 

13.40 

7 

10,500 

109 

86 

7.02 

1.94 

10,39 

9 

21,500 

156 

250 

9.  76 

2,85 

10.68 

10 

19,000 

120 

237 

lie  20 

2,18 

9.90 

11 

13,500 

257 

181.5 

13.95 

2„8? 

2.95 

12 

8,700 

228 

99 

13  c  20 

1.93 

5. 43 

1* 

13,300 

188 

195,8 

10.47 

1.35 

6.21 

15 

22,300 

129 

231 

8C  70 

2o56 

10,16 

18 

4,800 

125 

84o7 

4o50 

1.91 

6.05 

19 

760 

102 

66.1 

3.16 

1,93 

7.93 
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The  following  is  an  example  to  show  a  complete  calculation  of 
geomorphological  factors  studied  in  thi3  section: 
Watershed    S=l 

(1)  Area  A  »  100  sq<,  mi0 

(2)  Mean  relief 

cu: 


(i) 

/.6 


Figo  (7) 


R£igztve>    drea,    <*jfa 


A  «  100  sq»  mic 
H'  =  400  ft. 


Area  under  ths  hypsoasetric  curve 


-  J      7$  °W  =  °-^ 


(ii)  The  of  land  mass 

(ill)  The  mes     ::.sf 


r  2/6, 


(3)     Drainage  Density    D 
■  100  sq»  mio 

2E.  L  -  I066  mi- 

D  -  laity 


S  -   3,0  6o  mi -/so, mi. 
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(4)  Main  stream  slope   S 

Divide  the  main  stream  into  16  equal  sections ,  and  calculate  the 
slope  of  each  section.  Then,  a  weighted  main  stream  slope  can  be 
compxited.  The  calculation  Is  shown  in  the  folJor.ing  table. 


Table  VI 


10 


51 

1000  1000 
31,6 |31-6 


11 


12 


3 
1000  626 
31.6  25.0 


13 


4 
1250 
35,4 


14 


2 

2500 
50,0 


15 


16 


2 

2500 
50  0 


2 

2500 
50  0 


'  \i 


£(s?)     u  S/O.  I 


<?  •  /& 


n 


r 


n 


£(4**) 


.2 


g 


=  9.34   x  /o 


..<? 


(5)  The  shape  factor, 

A  -  100  3q„  mi 
Ac*0'  -  11,24  mi, 
A  0  -  29  >6  ado 


/ 


f- 


=  2.63 


113, 

VI.  MULTIPLE  CORRELATION 

Multiple  correlation  is  a  statistical  method  to  find  the  relatlc 
ship  among  one  dependent  variable  and  a  number  of  independent  variables . 
If  a  linear  relationship  exists,  then  the  method  of  fitting  is  to  make  the 
sura  of  the  squares  of  the  deviations  of  actual  observations  from  the  theore- 
tical linear  relation  to  be  a  minimum,  Thi3  is  so-called  the  method  of  least 
square „ 

General  Formula  for  Multiple  Correlation  (26) 

If     y.  «  i'th  observed  value  of  y 

y,  «  estimated  y. 
1  i 

-     <Z  *  htKti  +  &*Kai  +  b3X£  +  -  -  -   V-  ^  x*.- 
k    «  number  of  independent  variables 
n    •  number  of  observations 

a,  b^ii  bg,  b«j,  b, ,  ~-.--,«.c,„k    are  constants  obtained 
by  minimizing  the  sum  of  squares. 
Then  _.  />„         a^ 

Q=  £(Y<  -  K  ) 

*  £"(%•  -  a,  -  kX,;  ~  b^Xsi  -  63X31  -  •  •  *  -  bKK,ti) 

Taking  the  derivative  of  Q  with  respect  t®  a,  b. ,  b„9  b,,  -  -  —  b, , 

J-  <i  J  K 

and  equating  them  to  aero. 


ff  =    -  3(£?(,i  Vi  -ez,£x*i  -b.ZX*'  '  b*£x 


3  «% 


••*W b^Mx, 


> 


#&  -  -2<fefcc*/f -<aT**."   -  b,&*,iXm£  - £*£&)     s 
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Then,  we  havef 

an  +htgic,i  +6*€x*i  -/-  ■  •  m  +  k>*£Xx*'    -   £"K 

a£K,i  /  t>t£X*i    +  b££X,;X*  *•••*-  £*£*,*%/    ■   ^X*// 


The  first  equation  can  be  written  as 

^  ■  f  -b,K,  -  b£Xz  -63K5-   •  •  •  -  bK  X* 
substituting  into  other  equations,  we  obtains 


(7'6>X.~bJ*-  •  '  b*.XUL)£X,i+i>,£&?i  +  *&Klt-Xai  /•  .  •  •  +  J*s£j(l{  X*:    »  £X,t  yz 


(jt-kXt'biXt  ■•;  £«%.)£&/  +  t>XXaXii+bi£ x£i  +  •  *  •  +bk*teXKi    «  *X*iffd 


[Y -b^b£r£  .  <  -btx^Xti  +b,g*li;x,i  +M*£*t*X*i  *"■"  +6i£-X%    * «T<V  ^; 


By  rearranging,  these  equations  become: 


k£Xtfa-Z)+/>s£x,i{xtt"&)+'+6*£Xtifat>-j?M)  *£X*(tji  -p) 
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b&Kai  (x,i  -  X, )  y-  t^gX* (x2i  -  Xa )* "' *  b^E X^'Cx*;  -  X*  )    -  £ Xz<  ( yi  -  y) 


b,£XKi(Xti-X,)+ik£Xv:fa<-1?i)+  •  •'  +  b*£xtt£(x*i-x.k.}  *  &X*i  (<j)  -y) 


by  introducing  new  symbols 


S(x,z)  *£*%  ~-^a    -  £Xu(Xu-X) 

S(x,Xz     -  gXtiXzi  -  %£dg£*i    *  £Xd(X*i-  X* 

=  r  X,i  (x,i  -  %) 

* 

Then5  the  general  equations  of  multiple  correlations  are  as  follows: 

a,  »  jT  -  £, £  -£*Xe  -b3X~s~    ■    •   •   -  b^Y* 
bSfrf)+*kS(xtXg)  +6,S(x,X3)+'   ■  -+4*S<X,Xm)   *  SC**Y) 

btS(x,x*)  +  &3*Z)  +  ^S^e^  +  '  ■  •  -t-bkS(xzx*)  -  sc**y) 
btS(x3x,)  +  b*S(x3xe)  +  63S(if)  +  •   •  -^SfeaiJ  *  S(xsy) 

* 
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l       General  Formula  for  Estimating  0f  Variance 

Tha  estimate  of  tho  variance 

where  rf-(&  +  0        -  degree  of  freedom. 
The  sum  of  tha  square  of  the  deviations 

-*z(y  i  -  <2  -  b,X,i  -  £*A£«  -  •  •  •  -  htc  Xecf' 
-  £  iq;  ((ai  -co  -b,  X,;.  -  b&Xz;  -  b3 X&  -  »  •  •  -  hM X^O 

"~   cliffy  =c£-i>,  ^4  -  t&  Yet    -4^3* £>*:**<  ) 

-  b,  £X»l  (ifl-*.-b,  Xt ;  -  &9  Xzl  -b3  *W  - b*  XAc  ) 

■-bsc&titifa  -*>  -b.Xti-b***  -t>sX*i  -    •  •     "4e  X*<) 


All  terms p  except  the  first  one,  in  the  right  side  are  equal  to  aero 
since 

Then,  the  sua  of  square  of  the  deviations  can  be  expressed  as 


'(^s  ~  f<-)*  "  &<ji  (qi  -*>-  hX'i  -  *>*  **<'  -  *   *   '  -  *****) 
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Substitute  the  relation 


["nusB 


s    SVV)  -  i?,  SfX/fy)  -  &a  S(xzy)-  •  ■ '  -hcS(x*y. 


Henea^  the  general  formula  for  estimating  of  variance 


Co  Results  of Correlations ^between .Peak  Runoff  and  Gecmorphologlcal 

Factors  of  Sssall  Watersheds  in  Indiana 


The  following  Table  lH   shows  the  re  stilts  of  correlation  between 
the  25«years  instantaneous  peak  runoff  and  the  geoaorphological  factors 
of  3ns.ll  watersheds  in  Indiana  „ 


*4 
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Tabic  VII  -  Multiple  Correlation  between  Predicted  25-Years 
Instantaneous  Peak  Runoff  and  Geomorphologieal  Factors 


Runoff 


Geomc  rphclogical 
Factors 


A  f 

A  S 

A  D 

A  H 

ADS 

A  D  H 

A  D  f  E 

A  H  D  S  f 


Regression  Formula 


Q  -  30.29  A0"8471  f1"9245 

Q  -  0,0005716  A2  6792  S1*5829 

Q  »  47.13  A0'1994  D1  9847 

Q  -  0.03993  A0  7229  H1'7396 

Q  -  16,*5  A1*4179  D1"3292  S°°"S5 

Q  -  0,3432  ACo3925  H1"2732  D0"8601 

Q  -  0,8164  A0o3°36  H1"0936  D°°9409  S°  4098 

Q  -  Co05363  A0,9715  h°°7844  D°°8234  f°^i6°  s°  5903 


^Standard 
Derlatloc 


0*394 

0305 
0*296 
0  234 

0.2375 
0,213 
O0ZL8 
0,190 


*  In  Log-Unit 


It  is  quite  interesting  to  3ee  from  tha  table  that  the  larger  tha  number 
of  geomorphologieal  factors  put  into  multiple  correlation,  the  smaller  the 
standard  deviation,  That  is,  the  more  geomorphologieal  factors  we  considered, 
the  more  accurate  predicted  runoff  3&n  be  e:<pected6 

The  last  regression  formula  in  Table  VII 

Q  -  0,05363  A0"9715  H0"7844  D0'8234  f0"4160  S0*5901 
with  the  least  deviation  is  of  course  the  best  expression  obtained  from 
multiple  correlation,  If  all  the  geomorphologieal  factors  with  their 
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powers  together  are  combined  as  a  basin  characteristic,  B9 
Bx  -  a0"9715  H°o78^  d0c8234  f0„U60  sO,5901 

then,  the  regression  formula  can  be  stressed  as 

Q  -  C  BT"      C  »  constant 
which  is  the  equation  of  a  straight  line  in  log-log  paper.  Therefore, 
this  regression  formula  was  plot"     -a  a  theoretical  straight  line  was 
fitted  by  the  least  square  method. 

Figure  (8)  shows  the  result  of  plotting  the  basin  characteristic 
against  the  25~yoars  instantaneous  annual  peak  runoffs  It  seems  that 
the  deviation  of  actual  observation  to  the  theoretical  line  is  email, 
and  a  straight  line  can  be  fitted  quite  well.  A  95%  confidence  limit 
(26,  27)  *»s  calculated  and  plotted.  That  means  the  probability  of  any 
observation  falling  between  these  limits  is  95/S* 

Do  Construction  of  Correlation  Chart  for  Predicting  Future  Flood 

The  regression  formula  obtained  from  multiple  correlation  can  be 

used  to  calculate  the  flood  discharge  if  the  geomorphological  factors 

are  known^  and  also  tha  Figure  (8)  is  a  means  that  can  be  used  to  find 

the  flood  discharges  This  requires  knowing  the  different  power  of  each 

geomorphological  factor  before  the  flood  discharge  can  be  obtained-. 

There  is  a  somewhat  more  convenient  way  for  the  design  engineer 

to  convert  the  regression  foraula  to  a  chart  frora  whieh  the  predicting 

flood  discharge  can  be  read  directly  simply  by  entering  the  known  geo= 

morphological  factors  into  it„ 

The  way  to  convert  the  regression  formula  to  a  chart  is  as  follows: 

1)  Since  the  regression  formula  has  the  form 


Q=  efi'kB'c*D'iE 
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Figure  8 

Regression  Line  for  Theoretical  25  Xeara  Instantaneous  Feak 

Discharge  Against  Ea3in  Characteristics 
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95%    CONFIDENCE     INTERVAL. 

UCL=  UPPER    CONFIDENCE  LIMIT. 
LCL  =   LOWER     CONFIDENCE   LIMIT. 

X        * 

REGRESSION    LINE    FOR    THEORETICAL 

25    YEARS    INSTANTANEOUS     PEAK    DIS- 

CHARGE   AGAINST  BASIN    CHARACTERISTICS. 

"1,000  10,000 

BASIN  CHARACTERISTICS 

0.9715         0.7844  0.8234             0.4160 

(AH  D                       f 


0.5901 
S  ) 


A=DRAINAGE    AREA  (Sq  mi). 


D=  DRAINAGE    DENSITY  (  crib") 


MAIN     STREAM       SLOPE      (     S    X    10     ). 


Sq.mi ' 


H  =  MEAN     RELIEF  (Feet). 


f=  SHAPE      FACTOR. 


FIG.  8 
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2)  It  is  assumed  that 

Ht  *  A"3bCc   -    /f,Ca 
Ms  -A~e'c'i  Dd-  fia£>d 

Q,  &/?"■ Bb C"  Od £'   -     Krf3E* 

3)  The  above  equations  can  ba  plotted  when  Bt   C,  D„  Ep  are  used  as 
parameters «  They  are  shown  as  follows: 


A 


Fl«,  (9) 


tfi 


# 
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4)  The  above  figures  can  be  combined  as  shown  below: 

f   i- 


Fig  (10) 


As  shown  in  Figure  (10),  the  A&  is  in  arithmetic  scale.  This 
can  be  transformed  to  a  corresponding  Aa  scale  such  that  the 
area  A  can  be  introduced  directly  into  the  chart  and  join  Kith 
other  factors  to  find  the  flood  discharge » 


According  to  the  above  procedures  the  regression  foraula, 

Q  -  0,05363   A0  9715  H0'7844  D°°a234  f0**160  30-5901 
can  be  plotted  as  a  chart  for  the  convenient  use  of  the  design  engineer 
to  predict  the  25  years  instantaneous  peak  flood  discharge  for  highway 
drainage  structure  design  in  Indiana*.  There  are  two  charts  shown  in 
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Figure  (11)  and  Figure  (12).  Figure  (11)  ia  prepared  for  the  waters. 
less  than  100  square  mile3  and  Figure  (12)  is  prepared  for  the  watershed 
less  than  300  square  miles „ 

Two  examples  are  presented  as  to  show  how  to  use  the  charUto  obtain 
the  25  year  instantaneous  flood  discharge, 

a)  Watershed  No.  S-6 
Watershed  area    (A)     62,9  sq„  mi. 
Mean  relief      (H)     110   feet 
Drainage  density  (D)      S   mi<,/Gq0mi0 
Shape  factor     (f)      1.12 

Mean  stream  slops  (s)     21,10  x  10^ 

It  ia  shovm  as  the  dotted  line  in  Figure  (11)?  the  flood 
discharge  ia  read  from  the  chart  directly  as  38750  efs.p 
while  the  flood  predicted  by  the  frequency  study  is  3,300  cfa~ 

b)  Watershed  No*  S-=l 
Watershed  area  (A)  100  sq.  ado 
Mean  relief  (H)  216  feet 
Drainage  density  (D)       10.66  mio/sqomi,. 
Shape  factor  (f)       2o63 

Main  stream  slope  (S)        9o84  x  10 

This  is  shown  as  the  dotted  line  in  Figure  (12);  the  flood 
discharge  is  read  from  the  chart  directly  as  12,450  efs., 
while  the  flood  predicted  by  the  frequency  study  is  11,800  cfs.> 
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Figure  11 

25  Y©ars  Instantaneous  Annual  Peak  Runoff  for 

Watersheds  Less  Than  100  Square  Miles  in  Indiana 


FIG.  II 
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Figure  12 

25  Years  Instantaneous  Annual  Peak  Runoff  for 

Watersheds  Less  Than  300  Square  Miles  in  Indiana 
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VII.  CONCLUSIONS  AND  DISCUSSION 

lo  Table  IV  gives  the  predicted  flood  of  25,  50,  75,  100  years 
frequency  for  32  watersheds  in  Indiana, 

2,     Figure  (8)  is  a  plot  of  the  regression  formula 

Q  -  Oo05363  a°°9715  H°°7844  d°°8234  f0^60  s°°5901o 
This  can  be  used  to  find  the  25~year  peak  discharge  if  all  the 
geomorphological  factors  are  known  and  to  have  the  basin 
characteristic  calculated 0 

33,  Figure  (12.)   and  Figure  (12)  are  the  charts  used  to  find  the 
25=year  peak  discharge  for  the  watershed  less  than  100  square 
miles  asisd  300  square  sixes  respectively 0  From  the  two  examples 
shown  in  the  last  section,  it  is  learned  that  the  discharge 
read  from  the  chart  is  very  close  to  the  discharge  predicted  from 
statistical  frequency  analysis* 

4c  The  upper  and  lower  confidence  liiaits  in  Figure  (8)  mean  that 
the  probability  of  any  observation  falling  into  this  Interval 
is  95£*  It  ajay  be  safer  to  use  the  upper  limit  line  as  a  base 
for  design*  Then,  it  will  have  97*  5£  confidence j  that  iss  there 
is  only  2*5#  probability  that  the  flood  will  be  greater  than  the 
predicted  one  within  a  certain  frequency  of  25  years* 

5n  Therefore,  to  use  the  upper  limit  line  provides  safer  design. 
The  procedure  to  find  the  safe  design  discharge  is  as  follows: 

a.  Find  the  watershed  characteristics:  A,  H,  D,  f,  S* 

b.  From  Figure  (11)  or  Figure  (12),  read  the  25-years  instan- 
taneous annual  peak  discharge,  Q,  in  cfs* 
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Co  Go  back  to  Figure  (3)j  point  out  the  discharge  in  theoretical 
line,  and  go  straight  up  to  the  point  of  upper  limit;  then, 
read  the  corresponding  discharge „ 

60  This  report  contains  only  the  first  two  approaches  on  runoff  study, 
the  statistical  and  geomorphological  study.  It  i3  ba3ed  on  the 
assumption  that  the  rainfall  pattern  is  homogeneous  over  this 
area,  and  the  geoniorphological  factor  turns  out  to  be  the  only 
factor  affecting  the  amount  of  runoff 0 

?0  The  hydrograph  study  will  consider  the  different  rainfall  patterns 
and  their  relations  to  their  producing  runoff 0  The  theoretical 
instantaneous  unit  graph  and  th®  way  to  get  the  synthetic  hydro- 
graph  are  the  main  purposes  for  the  next  step  in.  the  study,, 


123, 


BIBLIOGRAPHY 

lo  Bow,R.  RoWnsoa,  "Hydrologic  Data  for -Highway  Design",  Amer.  Geoph      Loo 
Trans,  V.  28,  pp.  739=741,  Oct.  1947. 

2,  Geyer,  J.  C,  "New  Curve-fitting  Method  for  Analysis  of  Flood  K-cords"  Trans 

Amsr.  Geophys.  Union,  1940,  II,  pp.  660-663* 

3.  Gumbel,  E.  Jo,  "Probability  Interpretation  of  the  Observed  Return  Periods 

of  Floods?,  Amer.  Geophys.  Union  Trans.,  V.  21,  pt.  3,  pp.  336-350,  1941, 

kc     Gumbel,  B„  J.,  "On  the  Frequency  Distribution  of  Extreme  Values  in  Meteorological 
Data*,  Bui.  Amer0  Met.  Soc,  V.  23,  pp.  95-105,  Mar.  1942. 

5.  Gumbel,  E.  J.,  "On  the  Plotting  of  Flood  Discharges",  Amer.  Geophys.  Union 

Trans.,  ¥.  24,  Pt.  2,  pp.  699-719,  1943o 

6.  Gumbel,  S.  J.,  "Floods  Estimated  by  the  Probability  Method" ,  Eng.  New-Reo., 

¥  134,  pp.  97=101,  June  14,  1945 » 

7c  Gumbel,  Eo  J.,  "The  Statistical  Forecast  of  Floods",  Ohio  Water  Res.  Bd=.  21. 
pp.,  1943. 

3.  Gumbel,  E„  J.  "Statistics  of  Extremes",  New  York,  Columbia  University  Press. 
1953 . 

9o  Uc  So  Depto  of  Commerce,  "Probability  Tables  for  the  Analysis  of  Extreme- 
Value  Data",  National  Bureau  of  Standards,  Applied  Mathematics  Series  22B 
1953  o 

10o  Uo  S0  Dept.  of  ConsBerce,  "Statistical  Theory  of  Extreme  Values  and  Some 
Practical  Applications",  (A  series  of  lectures  by  Emil  J.  Gumbel), 
National  Bureau  of  Standards ,  Applied  Mathematics  Series  33,  1954. 

11.  Benson,  i-knuel  A.,  "Use  of  Historical  Data  in  Flood-Frequency  Analysis", 
Amer.  Geophy.  Union  Trans.,  V.  31,  pp.  419-424,  June  1950. 

12-.  Chow,  Yen  Te,  "A  General  Formula  for  Hydrologic  Frequency  Analysis",  Trans. 
Amer.  Geophys. Union,  Vol.  32,  No.  2,  April  1951,  pp.  231=237. 

13o  Whialer,  Benjamin  A.  and  Smith,  Charles,  "The  Estimation  of  the  Frequency  of 
Rare  Flocd3",  A.S.C.E.  Journal  of  Hydraulic  Division   Apr.  1957. 

Rows,  R.  R„,  Long,  Go  L»,  and  Royce,  To  C,  "Flood  Frequency  by  Regional 
Synthesis",  Trans.  Amer.  Geophys.  Union,  V.  33,  pp.  879=334,  1957. 

15  >  Langbein,  Wo  B.,  "Topographic  Characteristics  of  Drainage  Basin,  U.S.G.S. 
Water-Supply  Paper  963-C  155  pp.,  1947. 

16.  Strahler,  An  No,  "Hypsometric  (Area-Altitude)  Analysis  of  Erosions!  Topography 
Bui.  Geol.  Soc.  Amer.,  V.  63,  pp.  1117-1142,  1952. 


129. 


17.  Strahler,  Arthur  N.,  "Quantitative  Analysis  of  Watershed  Geomorphology'1 , 
Technical  Report  No.  13,  Dept,  of  Geol.,  Columbia  University,  1957, 

13.  Benson,  M.  A.,  ,:Channel-Slppe  Factor  in  Flood-Frequency  Analysis,  Jour. 
Bydr.  Biv.  A,S.C.E. ,  April  1959. 

190  Kinnison,  H.  B.,  ,:Flood-Flow  Formulas,  Jour.  Boston  Soc.  Civil  Engr.;  V.  33, 
pp.  1=19 s  Jan.  1946. 

20.  Yule,  Robert  B.,  "Bridge  v;aterway=Area  Formula  Developed  for  Indiana, " 
Civ.  Engo  V.  20  pp.  26,  73,  Oct.  1950. 

21 o  Bigwoed,  Bo  L.  and  Thomas,  M.P.A.,  "  A  Flood-Flow  Formula  for  Connecticut", 
U.S.G.S.  Circular  365,  16  pp.,  1955. 

22.  Dooge,  J.C.I.,  "Rational  Method  for  Estimating  Flood  Peaks",  Engineering 
y.  184*  No.  4774,  Sept.  6,  1957,  pp.  311-313,  No.  4776,  Sept*  20,  pp. 
374-377o 

23o  Gr^en,  Ao  Ro  and  Hoggstt,  R„  E„  "Floods  in  Indians  Magnitude  and  Frequency", 
U».  So  Dept.  o?   the  Interior  Geological  Survey,  I960. 

Dairysple,  Tate,  "Flood-Freouensy  Analyses",  U.SoG.S.  VIater=Supply  Paper 
1543-Ao 

25o  Taylor,  A.  B.  and  Sohwary,  H.  E.,  !SUnit~Kydrograph  Lag  and  Peak  Flow  Related 
to  Basin  Characteristics4/  Trans.  A.G.U.,  Vol.  33  p  pp.  235-246,  1952. 

26.  Bennett,  Co  A.  asd  Franklin,  N.  L.,  "Statistical  Analysis  in  Chemistry  and 

the  Chsadsal  Industry'^  Jehra  YJilsy  &  Sons,  Inc.  1954. 

27.  Burr,  I,  W„,  "Engineering  Statistics  and  Quality  Control",  McGraw-Hill  3ook 

Company,  Inc.,  1953* 


